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This Final Report is submitted in accordance with the require- 
ments of Contract No. NAS8-31629, '* Analytical Study of Inside-Out 
Gimbal D3mainics . ** 

This report consists of two volumes; 

Volume I - Analytical Study of Inside-Out/ 

Coincident Gimbal Dynamics 


Volume II - Appendix 
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INTRODUCTION 


Examination of candidate experiments projected for the joint 
European/ American Spacelab program indicates that a considerable 
number of these require pointing stabilities of + 1 arc-second 
or better which far exceeds the pointing stability capability of 
the Shuttle Orbiter. In addition, a number of experiments require 
accurate slewing, i.e., to within + 1 arc-second of a defined tra- 
jectory, in order to track a point on Earth with great precision. 

It is, therefore, clear that a pointing system is needed in order 
to meet Spacelab experiment requirements. To this end, the Euro- 
pean Space Agency (ESA) under the auspices of their ESTEC branch, 
working in conjunction with Dornier Systems, is developing an In- 
strument Pointing System (IPS) which is projected to meet Spacelab 
experiment requirements. The IPS being developed is a three axis 
gimballed system which mates with the back end of a prospective 
payload resulting in a large mass offset (i.e., distance between 
telescope C.M. and glmbal intersection point). The effects of the 
large mass offset are counteracted by a suspension system which 
isolates the IPS from disturbances such as man motion, RCS thruster 
firings, etc., that arise in the orbiter. The Bendix Corpor- 
ation, under contract to Marshall Spaceflight Center, has been 
conducting studies on the IPS for the past year. The objective 
of these stu d ies is to determine the performance capabilities of 
the IPS in order to aid NASA in their support of the ESA IPS design 
by maTn'ng the findings of the study available in a timely manner. 
Modifications that would result in overall improved IPS performance 
could be then Incorporated in the IPS design in a cost effective 
manner resulting in an experiment pointing system that will meet 
the broadest possible experiment requirements. 

This report describes the study results obtained to date on 
the performance capabilities and limitations of the IPS as present- 
ly conceived by ESTEC and offers some suggestions of design modifi- 
cations that if Implemented will result in overall improved IPS 
performance. Since the design and configuration of the IPS has 
been modified over the period of time this study has been performed, 
a portion of the study has been performed with the Inside-Out Gim- 
bal (lOG) configuration which was subsequently updated to the pre- 
sent Coincident Gimbal System configuration. However, due to the 
similarity of the two systems, the results obtained for the lOG 
also apply to the Coincident Gimbal System. 


This report consists of two volumes. The first, which is the 
main body of the report, describes the analyses and simulations per- 
formed during the course of the study and the conclusions and re- 
commendations resulting from these investigations. The second 
volume, which is an appendix to the main body of this report, pre- 
sents all of the eigenvalue data obtained upon which the results 
of the IPS system stability analysis presented in section 3.0 are 
based. In addition. Coincident Gimbal earth point tracking time 
histories that were not included in the main body of the report 
are presented in the appendix. 


2 . 


MODEL DESCRIPTION 


Two mathematical models describing two types of Instrument 
Pointing System (IPS) configurations were developed in the course 
of the study. The types of IPS for which mathematical models were 
developed are the Inside-Out Gimbal (lOG) and the Offset Gimbal 
System configurations. Only the lOG model was computer imple- 
mented since the Offset Gimbal System configuration was eliminated 
from contention as being a viable system configuration. 

The lOG model is a large angle nonlinear model consisting of 
three bodies. One body represents the shuttle orbiter and pallet, 
the second represents the lOG pedestal, and the third represents 
the instrument plus the inner (i.e., inertial) gimbal upon which 
the instrument mounts. Features of the model include; 

a. Full strapdown equations of motion describing the relative 
orientation between pedestal and instrument, 

b. Nonlinear Euler terms due to telescope rotation. 

c. Full six degree of freedom suspension dynamics between 
pallet and Coincident Gimbal pedestal. 

d. Complete representation of gimbal wire torques. 

e. Small and large angle tracking capability. 

f. Sensor and actuator dynamics. 

The Offset Gimbal math model is a large angle nonlinear model 
consisting of eight bodies. One body represents the shuttle orbi- 
ter and pallet, six bodies are used to represent the Offset Gimbal 
System configuration, and one body represents the instrument. 
Features of the mathematical formulation are; 

a. Offset gimbal flexibility characterization including flexi- 
ble interface between inertially stable gimbal and instrument. 

b. Full strapdown equations of motion. 

c. Nonlinear Euler terras due to various body rotations. 

f. Capability of mounting sensors on inertial gimbal or 
or instrument. 


In addition, a linear stability model was developed for the 
Coincident Gimbal by formulating the system matrix and obtain- 
ing the corresponding system eigenvalues. This model was used 
extensively to determine the stability boundaries of the coincl^ 
dent gimbal system configuration to control system and geometric 
parameLer variations • 

2.1 lOG Mathematical Model - The sections that follow will 
outline the derivation of the mathematical model used to deter- 
mine lOG pointing and slewing performance. 

2.1*1 Translational Equations of Motion - Referring to figure 
2“lj the following translational equations of motion can be bitten 
for the system depicted. 



(2-1) 


(2-2) 

^3e“V™3p3 

(2-3) 

3 



(2-4) 


j=l 

where ; 

F, Cj”l?2,3) = external forces applied to bodies 1 thru 
6 respectively 

^ci2 ~ compliance force between bodies 1 and 2, 

F “ hinge force between bodies 2 and 3 

m^(j=l,2,3) ~ mass of bodies 1 thru 3 

p.(j=l,2,3) ^ distance from origin oC arbitrary Inertial 
^ coordinate frame to center of mass of 

bodies 1 thru 3 
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= distance from origin of arbitrary inertial coordinate 


frame to the composite center of mass of the system shown 
in figure 2-1. 


Additionally, from the geometry shown in figure 2'~1 we have 




n =p 4-]^ 

*^0 3 


where! 


R. (j=l,2,3) = distance from composite center of mass to the 


center of mass of bodies 1 thru 3 


From the definition of the composite system center of mass 
Also from geometrical considerations 


^2 ’^ l ’’'^ 120‘'‘^12 


^3 V ^ 120 ‘‘’^ 12 '*'^ l '*^2 


where: 

R- = inertially fixed vector directed from the CM of body 
1 to body 2 when the system is in an unstressed state 

^12 ° relative linear displacement between bodies 1 and 2 

measured with respect to R-j^20 

r^ = distance from CM of body 2 to gimbal hinge point 

r^ = distance from gimbal hinge point to CM of body 3 

Substituting equations (2-9) and (2-10) into equation (2-8) 
and solving for yields 


Substituting equations (2-7) , (2-11) , and (2-4) into equation 
(2-3) and solving for the hinge force yields 


( 2 - 12 ) 


The compliance force between bodies 1 and 2 of the k lso“ 
lator can be written as 


^cl2k"^S2*‘^12‘*'^12 




where : 


{iJJ‘(k=l,...,4) 

3 ^ 1 o(k=l,..., 4 ) 

3 ^o(k=l,..., 4 ) 


distance from CM of body 1 to the lumped 
spring damper between bodies 1 and 2 on the 

til 

body 1 side of the k spring damper system. 
This vector is fixed in body 1. 


t h 

distance from CM of body 2 to the k lumped 
spring damper between bodies 1 and 2 on the 

t h 

body 2 side of the k spring damper system. 
This vector is fixed in body 2* 

12 

value of when system is unstressed. 


value of 



when system is unstressed. 


^12 ” spring constant between bodies 1 and 2. 

= damping constant of spring damper between 
bodies 1 and 2. 

The damping and spring constants are defined as diagonal 
matrices in the following manner 


kj2=‘‘li^g 

I>12=dlag 


fl2»"‘l2yf'l2^ 

pi2x’“l2y*“l2^ 


(2-14) 

(2-13) 
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For small angular rotations of bodies 1 and 2 the following 
relations apply 


In addition 



(2-16) 


(2-17) 

,, „12 
®kl“V\lO 

(2-18) 

fil2 ,, fil2 

*^k2 “2*®k20 

(2-19) 

kl0"^120'^^20‘^k20 

(2-20) 

8^^ — R 

^k20 ‘^E20"*k20 

(2-21) 


where: 


0^ = angular rotation of body 1 

©2 = angular rotation of body 2 

(0^ = angular rate of body 1 

- angular rate of body 2 


*^20 


= inertial vector equal to the distance from CM 
of body 2 to center of elasticity of spring 
damper system between bodies 1 and 2 when the 
system is unstressed 




inertial vector equal to the distance from the 
center of elasticity of the spring damper sys- 
tem between bodies 1 and 2 to the spring 
damper of that system. 
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Substituting equations C2~16) thru (2-21) into equation (2-13) 
and rearranging terms yields 


F 


[^12’^ ^‘^2"“i ^ ^E20'*'^“2''“l ^ ’^“lc20''‘^1^12 o"j 


“l2 


( 2 - 22 ) 


However 


F 


cl2 Z-j cl2k 
k=l 


(2-23) 


Substituting equation (2-22) into equation (2-23) and performing 
the indicated summations yields 


12 1 

^cl2"^12* r^l2'^^^®2'‘®l^^^E20'"^®1^120‘*‘^®2'‘®l^^ Xi ^^20 

^ k=l -* 

4 _ 

k=l 

However, by the definition of center of elasticity 

4 


Therefore 


E 


k=l 


\l2=^^12 


+4D 


12 


* [^12'*'^®2~®l^^*^E20"‘^i^l2(^ 

‘ j^l2‘**^‘^2"‘^l^’'^20““l^*^12<^ 


(2-24) 


(2-25) 


(2-26) 


V ? 

i { 


0 

Q 


0 

a 

0 

0 
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A force equation will be required in order to solve for e^ 2 * This 

equation is obtained by substituting equations (2-26) » (2-11) and 
(2-5) into equation (2-1) which results in 


M 


m. 


^e- 


^^(F^g+Fje)" 


-P 


(m2'hn3)ej2+4D-ej2+‘K^2' 


=12] 


+4D 


12 


|^(a)j^-u) 


2^**^E20'^1*^120 


0 ■‘■^^12* 


® 2 ^^^ E 20 ’*’® 1^^120 


} 


iT 


(V-2> 


(2-27) 


2.1.2 Rotational Equations of Motion - The following para- 
graphs will develop the rotational equations of motion of the 
system depicted in figure 2-1. In this development it will be 
assumed that the coordinate frames of bodies 1 and 2 are aligned. 
However, the coordinate frame of body 3 can be at an arbitrary 
orientation with respect to bodies 1 and 2. 


Body 1 - The rotational equation of motion for body 1 can be 
written as 

4 


T, +R, ,xF, + 
le 11 le 


k=l 


kl^^cl2k dt^'^l’^^l^ 


Substituting equations (2-16), (2-20) and (2-22) into equation (2-28) 
and doing the indicated summation eliminating second order terms 
yields 


T, +R„xF, = ^(J *10, )+d, ,, 
le 11 le dt^ 1 r 12 


(Wj^-W 2 )+kj^..* ( 0 ^- 62 ) 


■^<“l20+“E2o)=‘j“l2*^V“2>*®E20'^l’'‘'l2o] 

*12- [«^-®2>’®E20-"®l=®12 J-«\20™E20^*[^2-^12«12-^ <2-29) 
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Body 2 - The rotational equation of motion for body 2 can be 
written as 

4 




- I C 


■^12k'^l“V d7^-'2‘“2> 


(2-30) 


k=l 


Substituting equations (2-12), (2-17), (2-21), and (2-22) into 
equation (2-30) and eliminating second order terms gives 


■'HT*''2e«22’‘^2c* 


(m^-rt.2) 2 ">3 


M 




+k 


12 


m. 


(02-Oi)+4Rg2o*|^12* |^‘‘’2'‘^l^^S20~**’l^*H2^‘*'‘^12’ |^^2"®1^^‘^E20“^1**42^^ 

,m« r , ^ m (in +in„) „ 

IT ’^1’'^12'^^^E20^|^12*^X2'*^12*^1^'*' M **1 


(2-31) 


where : 


- external force acting on body 3 written in body 2 coor- 
dinates 

2 

= hinge torque acting on bodies 2 and 3 written in body 2 
coordinates 


= wire torques applied to body 2 written in body 2 coordinates 
w i 


0 

n 


where 

k^2=<iiag{4K^2z^“l20y^ ’ ^®l2z ^“l20x^ ’^^^2y^°^120x^ "‘■^12x^“l20y^ 

di2=diag{4ni2^(ai20y) “l20x^ '’*°12x^“l20y^ 

Tj^(j=l,2,3)- External torque acting on body 


I ; 


r^ - distance from hinge point to CH of body 3 written in body 
2 coordinates 
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Body 3 - The rotational equation of motion for body 3 can be 
expressed as 




(2-32) 


Substituting equation (144) Into equation (154) yields 




(m^^-hn^) 




M 




„ . V3 y3 . 

“V M ’^2^^12*^ M 


r2x(r^+r2) 


(2-33) 


where: 

Q 

^ wire torques applied to body 3 written in body 3 
coordinates 

jFg = external force acting on bodies 1 and 2, respectively, 
written in body 3 coordinates 

(a)^ = angular rate of body 3 

= distance from CM of body 2 to hinge point, written 
in body 3 coordinates 

Rewriting equations (2-27), (2-29), (2-31), and (2-33) in matrix 
form yields 


'^le+^ll‘''le='Jl*^1^12* <‘*‘l”^2>-^’^12* <®r®2> 
"^^^120‘*‘^20^ ’{^12‘(^^E20* ^*^l"“2^'*'^120*“^ 

■‘■^12* (fE20* ^®1"®2^‘‘‘^120*®:^J 

“^^^120'*‘®E20^ * 012*^12'^^12‘‘^1^ 


(2-34) 
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and 2^2 is the transfonmation from body 3 to body 2 coordinates and 
is equal to 


T = 
2 3 


^11 

^12 

®13 

^21 

^22 

^23 

^31 

^32 

^33 


m -1 

®ii 

^21 

^31 

T = T - 

2 3 3^2 

^12 

^22 

^32 


^13 

^23 



T =- 
H 


u) 


hT 


s*Hx} 


nT 


CO 


ns 


2 2 
s +2? CO s-Ko 

s ns ns 




■a)*+Kp 


a* ®3 
'®3+^I s 


J 


and 


“3=^“3“^c> 

e ;=( G 3 - 03 ^) 


(2-41) 


(2-42) 


An assumption inherent in equations (2-34) thru (2-37) is that 
the nonlinear Euler, centrifugal and Coriolis terras are negligible 
for bodies 1 and 2. Equations (2-34) thru (2-37) form the complete 
set of equations of motion of the system shown in figure 2-1, 

It now becomes necessary to define the control torques and 
the computation scheme for updating 2^3* control torques are 

given by 


(2-43) 

(2-44) 

(2-45) 


63^,= / “^.dt 


3c 


(2-46) 


where: 


K^.dlag[K^ 

• S - 'S 


(2-47) 

(2-48) 
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(2-49) 

% 


t, - rate gyro damping ratio 

2.1.3 Computational Scheme for Updating “ In order to 

complete the equations of motion given in equations (2-34) thru 
(2^37) , a method for computationally updating the transformation 
from body 3 to body 2 coordinates must be specified. Since the 
effects of girabal wire torques is to be evaluated in the course 
of the study it will be convenient to express ^T 2 in terms of modi- 
fied Euler angles and update the transformation by integration of 
Euler rate equations. In this manner, the Euler angles and rates 
will always be known, x^hich will be convenient for computing wire 
torques and gimbal friction x<h,ich depend functionally on Euler 
angles and rates. In addition, none of the gimbal angles approach 
77/2 rad (90 degrees) which eliminates the singularity problem 
encountered xfhen Euler angles are used as the strapdown computation 
algorithm. 

In order to properly compute the transformation between 
bodies 2 and 3 coordinates the rotation sequence between the bodies 
must be specified. For the IBS that rotation sequence is z, y, x 
when going from body 2 to body 3. In addition, since bodies 2 and 
3 rotate independently of each other relative Euler angles and rates 
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which are a function of these rotations must be computed in order 
to properly specify the transformation. Equations (2-50) thru 
(2-52) specify the transformation matrices going from body 2 to 
body 3. 


COS 0^ 

sin 0= 

--stn 02 

COS 0H 

O 

0 


1 2 . 


(2-50) 


0 cosfx. sin% yo 

0 -sin% cosfx, Bo 


C2-51) 


0 -stn 


C05 ry 

O 


sin 0 


C2-52) 


where 



= x,y. 

^o»^o*2o 

= x,y. 

X ,Y ,Z 

=■ x,y. 


= x,y, 


= x,y, and z axis of the pedestal fixed coordinate frame 

= x,y, and z axes of the outer gimbal 

= x,yj and a axes of the middle gimbal 

= x,y, and z axes of the inner or inertial gimbal 

= relative rotation of outer gimbal about the z axis 
= relative rotation of the middle gimbal about the x axis 
= relative rotation of the inner gimbal about the y axis 
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rr 

f 


However 


9^ - ©2.S 

n = %y - %y 


> 


F 


C2-53) - 


where 


0^^(k=2j3) = Euler angle about the z axis due to the k body 
inertial rotation 

<j)^(k“2y3) = Euler angle about the x axis due to the k^^ body 
inertial rotation 

th 

Clt”2,3) = Euler angle about the y axis due to the k body 
^ inertial rotation 

The telescope (i.e., body 3) inertial rate written in terms 
of Euler rates are given by 






-^3?. 



cos '^Y ^ 


■ si n % 


0 1 O 

Expanding equation (2*-54) yields 


0 


0 

cos 

% 

sin (Px 

-sin 


cos 



‘ 9 


' 0 ■ 


0 


K 




_ 0 


(2-54) 


V.x’ 


-cos^sinHy cosily 0 



fJ:,y 


sinfp^ O 1 





C05 cas^Y S/n 0 


'P 

L^SyJ 


(2-55) 


Solving equation (2-55) for the Euler rates gives 


© 3 / 


-sec^^s//»^^ 0 sec%'cos'^y 


- 3 ; 

^3X 

= 

cos V^y d SirtV'y 


^3y 

1 


_ "tein s/o 1 “ion Cos '' 


W 35 


The pedestal (i.e., body 2) inertial rates written in terms 
of Euler rates are given by 



COS 02 

“ ?i./7 ©2 

0 


1 

0 

0 


«*«• 

X 

1 


” 0 ' 


s/n 

cos 

0 


0 

cos % 





0 


0 

0 

i 


0 

sin ^ 

COS "Pk 


_0 

i 



Expanding and simplifying equation (2-57) gives 


1 


0 Cos 0^ -sinOg CosiPjf 


®2a 


= 

0 sin 6^ cosGgcos^j^ 


L 

ftt 

1 


1 0 Sio 




Solving equation (2-58) for the Euler rates due to body 2 


gxves 




2S 








02 ^ ^0 5 0£ iart VJc 1 


sin ©2 

COS 03 




a). 




U>: 




(aJ 


2Z 


( 2 - 


Cos ©2 sec iPx 0 



1 


M 


The relative Euler rates are obtained by differentiating 
equation (2-53) and substituting equations (2-56) and (2-59) into 
the results which yields 


- 

t - 
% ^ 


^Y — (u)^^ Cori ^ £*-^y Sfr> 0^) 

Wjx i<an % C05^Y “ f y 5CC Co 5 0^ - 0)^^ ?(f C J/r) 




The transformation from body 3 to body 2 is given by 


aT, 


COS 0 j — S ( o O 

Sir ^ Q~ o 

0 'T 1 


10 o 
O cos % - sin 

0 sin fy, f /I T 


[tOsV^y 0 


0 10 

-SmVy O C o s 


{2-U) 


Expanding equation (2-61) yields 



CO S SjCO.-T Vy “ 


* sm ©2 Cns 

C.OS *9y ■" 0^ "rjs V^y 


T ^ 
^3 

tin 05 COsV'j, 

COZ r./nS^y 

Con. Cos 

rin9f s 


h-a) 


- Co.S ' 

?y S,o 'Py 

Stn 'P^ 


cos 9^ rnsPy j 

- 


2.1.4 

Specification of 

Gimbal Wire Torques 

- The gimbal wire 

— 


torque characteristic is shown in figure 2-2, where 

^0z ” torque due to relative Euler rotation about the z axis (n-m) 

= Wire torque due to relative Euler rotation about the x axis (n-m) 

= Wire torque due to relative Euler rotation about the y axis (n-ra) 
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= Wire torque slope 

H^T " JuniP ill wire torque characteristic when relative Euler rate 
changes sign (n-m) 


Using the transformations given in equations (2-50) thru 
(2-32) the wire torques acting on body 3 (i.e, , the telescope) 
can be written as 




COS Vy 

0 



1 

0 

0 


- J 

1 


0 

'^WT3y 

- 

o 

i 

0 


0 

cos % 

sift 


O 

-h 


JwTSi: J 


sin 

0 

cos fy 


0 


cos 


n 

1 


0 

. _ 


(2-i3) 


Expanding equation (2-63) and rearranging terms yields 


T 


^ — 

-COS%Sini^y COS^y O 



^WT3y 

- 

Sinf^ 0 1 



T 


cos Cos sin*^y 0 




(2-64) 


Similarly the wire torques acting on body 2, i.e,, the lOG 
pedestal are given by 


T 


0 cos - sin 9 CJ>S ‘Px 



^T2y 

- 

0 Sin BjL COS0g CJDsfyt 



T 

'vyr2z 


1 

O 




(2-65*) 


Equation (2-64) and (2-65) combined with equations (2-34), (2-35), 
(2-36) , (2-37) , (2-43) , (2-60) , and (2-62) form the complete mathe- 
matical model of the 106. 


Table 2-1 lists a baseline set of parameters for the lOG 


4 

10^^=8.054x10 kg 

“l2x°“l2y°'’l22=59.93 n-sec/m 

102=1.95x10^ kg 

^12x=V^122-1“ 

012=2.683x10^ kg 

dj^2x“k^*96 (n-m-sec)/rad 

J^^=l. 551x10^ kg-tn^ 

d, „ =14.96 (n-m-sec) /rad 
12y 

fi 2 

Jj^y«6. 953x10° kg-m 

d-„ =29.92 (n-m-sec) /rad 

X^Z 

J^^=7. 543x10^ kg-m^ 

=31.25 n-m/rad 
12x 

J2x= 50 kg-m^ 

k. _ =31.25 ii“in/rad 
12y 

J2y=50 kg-m^ 

ki22=62.5 n-m/rad 

J2^=50 kg-m^ 

. 214x10^ (n-m-sec) /rad 

J, =2.903x10^ kg-ni^ 
3x 

Kj^^^=l. 972x10^ (n-m-sec) /rad 

J^^=2. 648x10^ kg“ia^ 

47x10^ (n-m-sec) /rad 

605x10^ kg-m^ 

Kp ^^=7 . 856x10^ n-m/rad 

R =15 1 +4.7x10'’^ 1 m 
11 X z 

4 

Kp^y=7. 00x10 n-m/rad 

R _^-0.375 1 m 

K =8.764x10^ n-m/rad 

E20 z 

P3z 

r-=0,375 1 m | 

J- z 1 

Ki 3^=1. 24x10^ n-m/sec 

r_=1.689 1 m 
2 z 

Kj- 3 y=l . 105x10^ n-m/sec 

j\ y\ 

5 

Rion=-^*72 1 -1.579 1 m 

120 X z 

Kjg^=l . 384x10^ n-m/sec 

W^.ji=125.7 rad/sec 

W =163.4 rad/sec 
ns 

?^=0.5 





) 


2,2 Offset Gjbnbal Mathematical Model - The sections that 
follow will outline the derivation of the mathematical model for 
the offset gimbal system configuration. Figure 2-3 is a schematic 
representation of the offset gimbal for which the mathematical model 
described below is derived, 

2.2.1 Translational Equations of Motion - Referring to fig- 
ure 2-3, the following translational equations of motion apply 




“ ^CIZ ^H23 ~ 

( 2 - 67 ) 

^134 " ^ 3/^3 

( 2 - 68 ) 


( 2 - 69 ) 


( 7 - 70 ) 


(2-7/) 

^7e " ^67 ^78 ' ^7/®7 

(2-72) 

^Se. ~ ^Z78 ~ 

(2-73) 

8 

if “ 'Vo 

(2-74) 

6 


(2-75-) 




where 

F. Cj“l»»*-»8) = External forces acting on bodies 1 thru 8 

j ® 

p (j=l, ,,,,8) = Location of CM of bodies 1 thru 8 with respect 

^ to some arbitrary inertial reference frame 

m. Cj-1, , . , ,8) “ Mass of bodies 1 thru 8 
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r 


“ Suspension force between bodies 1 and 2 
= Hinge forces between bodies 2&3, 4&5j and 6&7 

riXo flAl? Hd j . m. f 

respectxvely 

p = Location of system composite CM 

o 

F__ . ;F__^;F^_q = Interface forces acting between bodies 3&4, 

XDO J . / o c ^ j c o 

5&6, and 7&8 

From the geometry shown in figure (2—3) the following relationships 
apply 


ft 

= A " 


ft 

= A * 


ft^ 

- A ^ 

^3 

A 

-ft. " 


A = 



A ^ 

^ Po * 


A = 

-ft. ^ 


A = 

^ A " 

Rs 


^2 * * ^>2. 

^3 = 1^170 * ->7 

R-20 " ^ n * 

R 5 = ^ £.7 * n ^ * n, * r-^ + 

- R. * Rn:> ^-17 * f, * n, 

= R, > Ran f (- a ^ r, # r, ^ V c, ^ ^ r, , r-,^ 

Rg ^ /?, - .^,.-0 4 ■) r,^ 4 i ^ ^ r> -t * Ho * n, * r,^ 


(2-76) 

( 2 - 77 ) 

(2-7S) 

(2-79) 

( 2 - so) 
( 2 - 8 /) 
(2-SZ) 
(2-83) 

(2-S>4) 

( 2 -as-) 

( 2 - 84 ) 

(2-S7) 

(2-ffS) 

(2-89) 

(2-90) 
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f 


where 

^120 ” vector corresponding to the location of body with 

respect to body at the start of the problem with the isolators in 
an undeformed state. 

From the deflaition of the center of mass the following applies 


ft-’') 

j-' 

which results in 

ft-n) 

Substituting equations (2—83) thru (2—90) into equation (2-92), rearrang- 
ing terms and solving for gives 

■* ^ + m ^ 4- €|q^) + ( mj ^ * my 4- mg)fr| f Tjj^) 

^ ^ + jTij. + 4^ my ^ + Hij } 4- f ^ ^ * 'aJ 

+ (rnj4 m 7 + mj)rr 7 +rg) + (?«, + WpXr, r,^) * rn ^ r„ * r,^)] (2-'?3) 


Substituting equations (2-84) thru (2-90), and (2-93 into equations 
(2-76) thru (2-83) differentiating and using the relationship 




(a-94) 


results in 


^ I p ^ fl . 

/? ~ '"j 

e 8 


2-21 


0, " TT f'’ ~ ^ ^ ^3* I 


■ r^ ^ ^)j <• {r; ^ r'g} / .i rr>-j i- ( 


Pi ^ -pf ■^ Tr ^‘3 ^ ~TT' ^*^ *’'*^ 


r’ + ru,')(rri7* my") <- (r„ < r,j)mg ] 


[ fr; . r^ / 


+ (r^ ♦ ^4 ^2 + (^7 ■•■ 4-’,. » ; *■ ^ - 


I ^C~> - nt, „ m,^■m■2. , .. 

« <^ 7 ' ' 


^r, ■* >i) *7^(0* *'4' 2 wj 


® \ n 

■ f^[(rV+ '"j + (r'7 rj)^ /r>^ + (r^ ^■ r^^y.m^-i m^) *■ (r„ * r,^') J 




/O5 = 2 A f ♦ '^h'>(n * ^2) <''^3 ^ ^ 4 ^ ^'^’j + < '"r * 'a > ^ '"j 

"7f [(rV *- f's'. 2 tTr’i <• ^ r^fr'^ ^ ^ ^ 2 ) -^aJ 


= 7^2 li - 1 ‘ ^'"'' ‘ 'V + + 

^Ciy^r^'jAjf^l - 77 L + *• ^ + r, 2 .^ 

>' 

I ® t r .... V* .. .. ^ 

= 77- 2 Fje + Tfi 'n.c.j + (m, vir^^fr,' t n.' +■ (r^ + q) ^ mj i- ^ 

.. .. s ^ .. .. "I *”a / '• •• \ 

■*■ (i* 7 + m, ^ {^•* r",.,! ^1 frij j tf'ii'*' Hi ) 

>' j=< 

^ = ;;7 2 fj,. ■(- 77 ['^,e,^ + (m.^frij'ur, + Fj) ■« <r, <-^') 2 '”t ^ 

^ -1 
+ (r, f n; .'• . ^ 7 - +(r,frJZ-.m, f r' r,.J Zj rr.- i 


^■^v 



POOR QUALITY 


Substituting equation (2-102) into equation (2-73) and rearranging terms 
gives 


I 7 y 3 

^ -f- + «i) + ( »^3 + ^ § "’j 

J“* J=# * J 


** 5* 6 7 

•*'(rj,+ f;) + (7; + TglS (r^ + Oa) £ "’j + ^ + Hi) S m; 1 (2-/03) 


j=i 






Substituting equations (2-101) and (2-103) into equation (2-72) yields 
6 ^ 


^47 “ - (»«7 + m,) £ F, ] - ^/2. +fm,,^«7,)(A; + r;) 

j-' • 1=/ 

^ ^ ^ 

+ + ^ rr.) (i\ ♦ r)o) ^ j" TT (^ii +^71) ^ 


*(r, 


Substituting equations (2-104) and (2-100) into equation (2-71) yields 

a . _ 3. 


^£:4'M r SfnjS F,'e “ £ »”i£ ^ ( "’i ^ 

^ t i jrj J -* J=6 J=I '' j = t J-^ 

Smj (r; + ^ )S mjb ' + fUg ( rj, + r^) j] 

j - I J- 1 J" * 


Substituting equations (2-99) and (2-105) Into equation (2-70) gives 


S' a a A 


f^45=H “m if ^*11 +fm, + mjl(r; + r;)4-(ri4rL,)S'^3 

4 4 P 

fj.) -^- { )[{r7 + l'-^) mj <- (m? + nOg)(r^ -f r'^) + mg(r), -» r,^)]| 

1-1 >•■ j=6 


Substituting equations (2-98) and (2-106) into equation (2-69) yields 
3-S. e3_.^,_e_._ _|. 


=M C£wj^^'e “ £ ) +(^ + Itij] 

' ' j=t j=4 jH i=/ " *• j-4 J=l 

+ (S mj)r(rj+r ^)2 m 3 + (ry + rg)Smj + (m7 + mg)ff^ t Ho) ■^ rTtgCr,, *• f 7 ,)l} (2-^07) 

}~ * j-5* i“6 


Substituting equations (2-97) and (2-107) into equation (2-68) results in 

S 8 _ * r ^ 


%3~ + " (f5i + F;2)Smj] ~]^{(2mj)[»n,e,^ + fm, + rr>,)fr;v.rj)l 

j-3 ' ■ j:3 J" ■ 

^ S 6 ^ 

^ Crr?,+ rn^ i(r3 + "jZ m3 + (<v + q}2 Wj ^ fr; * r£ rr>j + (rn,-# mg)(r.j + rjg) -* fOg fr^, + )]j (Z-/08 




The Suspension Force derived in exactly the same manner ^ 

as described in section 2.1.1 and is given in equation (2-26) which 
is repeated here for convenience 

fcf2 ^ 4D-r(wj-f^,)y.Rp2j-qD*(uj,xR,.„) +^K-e,„ fo, \ (2-26) 

2o2.2 Rotational Equations of Motion - Referring to figure 
Z-^3y the follo^^ing rotational equations of motion apply for the 
8 bodies depicted 







^ if 

- ji 

■k 


) 






(2-m) 





*'%3WT ^ 

- Sa' 

^ p 

9 ^ Hliik 

'I 

r M^i 

^ 0/23 ”3 Jah 

_ d 


j (2-//0) 


' 



+ Rjj . 


-!• r- X /%23 ■'a 



cl 

dt 

^^3 

' oJj ) ->- dj^ 



^34 '( dj- 

©^) “ 7^2 

(2-///) 






f n,>i%4 +•> 

^ ^4C ~ ‘ 


d 

?T 

(x,. 


<(uj^- 

.w,) + 



(2-//2) ; 



~^5t/yT ^ 


’ flxF^^o+r- ^ 

^ (» *“ ■ 


d 

dt 

i^s 

» CJf) + dsiS" 

iui^. 




(2-//i?) 

1^. 


+ %7wr 

+ Rf,6 



*RiiA7 ~ V 

K1 ’ 

d 

d^ 

(^A- 

•6>J < d^. 

(iOn‘ 




(2-//4'l 


"^67/ 

" "^7wr 

+ Rfj 



'-Frrs^: 

h 

d 

d-* 

(Tr 

. /■.«),,) 4- dyg. 

■(to. 

Wg) 

^ l<7g *(07 

-0fl)-Tyy 

l7-i/y)\ 



^8. 

^8<5 

’<^SPc 

N2 ^ 'X7|E? 

d , 

■Ug) + cJyg * ('Jjg: 

CUy) + 

'(dg- 

0,) 


(2-R6) 



where 

12 

®kl- *'cX2fc’ ^H23! '‘l34’ ^HA5’ ^67 F„g ara defined In sections 
2,2.1 and 2.1.1 

R. . Cj“l»...8) = Distance from CM to the applied external force for 
the jth body 

T„ -jT, _ = Frictional torques between bodies 2&3, 4&5, and 6&7 

23f 45f 67f respectively 

J. (j=l, . . . ,8) = Inertia tensor of jth body 
3 
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respectively 


C I 

3^2 1 2 Constraint torque applied by body 2 on body 3 written in 
body 2 coordinates 

I - = Constraint torque applied by body 2 on body 3 written in 
body 3 coordinates 

_M,|c = Constraint uorque applied by body 4 on body 5 written in 
body 4 coordinates 

cM^L = Constraint torque applied by body 4 on body 5 written in 
body 5 coordinates 

= Constraint torque applied by body 6 on body 7 witten in 
body 6 coordinates 

= Constraint torque applied by body 6 on body 7 written in 
body 7 coordinates 

T^^,Tg^,T = Gimbal torques applied by the gimbal torque motors 

p* Ty about the a , P , and Y axes respectively 

Z IX- ^ 


In order to obtain a complete set of equations of motion for 
the offset gimbal system a translational equation of motion is re- 
quired to be added to equations (2-109) thru (2-116) . Equation 
(2-67) will be used to fulfill this requirement and is repeated 
here for convenience 


^23’‘^C12^^H23 “ ^*^2 (2-67) 

Equations (2-68) and (2-109) thru (2-116) form a complete 
set of equations of motion for the Offset Gimbal system configura- 
tion. 


Substituting equations (2-13) , (2-16) thru (2-21) , (2-26) , and 
(2-103) thru (2-108) doing the indicated inertial time differentia- 
tions, summations, simplifying, rearranging terms, and putting the 
equations in matrix format yields the following set of system equa- 
tions of motion, all of which use the following convention for the 
mass stmunations: 





( 2 - 117 ) 
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+ 4 . ) i- k* -90 ”‘^f^?, 7 (j+ Re 5 o)‘(^'^^£ 2 o*''^i 
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where 


= Constraint torque applied by body 2 on body 3 about the 

2 X and y axes written in body 2 coordinates 

2 M 2 I = = Constraint torque applied by body 2 on body 3 about the 

3 X and y axes written on body 3 coordinates 


c I ^ 

= Constraint torque applied by body 4 on body 5 about the 
4 y and z axes written in body 4 coordinates 


5^4 


y,z 


Constraint torque applied by body 4 on body 5 about the 
y and z axes written in body 5 coordinates 


7«6l 


x,z 


Constraint torque applied by body 6 on body 7 about the 
X and z axes written in body 6 coordinates 




x.z 


Constraint torque applied by body 6 on body 7 about the 
x and z axes written in body 7 coordinates 


The definition of the various matrices indicated in equations 
(2-118) thru (2-126) is given at the end of this section. 


The reason the constraint torques only appear about two axes 
reflects that each of the gimbals have a single rotational degree 
of freedom. In addition, the girabal friction and wire torques 
only act about the cixis of the respective gimbal rotational degree 
of freedom as indicated in the above equations. 

Inserting the gimbal constraints, the following relationships 
result: 





Equations (2-146) thru (2-147) are used to compute the rela- 
tive orientations of the various bodies or they form the strapdown 
equations of motion. 

In order to complete the equation of motion development for 
the offset gimbal system configuration, the constraint torques 
must be eliminated from equations (2-119) thru (2-124). 

The constraint torques may be eliminated by solving equa- 
tions (2-120) , (2-122) , and (2-124) for these torques which only 
act about two axes. Hence the constraint torques are given by 
the expressions that follow; 
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(2 -ISO) 


It should be noted that in the derivation of the gimbal con- 
straint torques only the linear terms were considered, and every- 
thing that follows from this point on considers only the linear 
terms of the equations of motion. 

Substitution of equations (2-148) , (2-149) , and (2-150) into 
equations (2-119), (2-121), and (2-123) respectively results in 
a set of system equations that can be written in the following 
manner : 
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2.2.3 Derivation of Glmbal Torques - In order to complete the 
mathematical model for the Offset Gimbal, the gimbal torque equations 
must be specified. Using the transformations shown in equations 
(2-136), (2-13fl) and (2-140) the desired control torque can be written 
in terms of the gimbal torque in the following manner: 


cosJTy O ' sin^yCosyS,. I 


Tc = I 0 i 


f! ^y 0 cosS‘y cos^y J LXt*':. 
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Solving equation (2-350) for the gimbal torques gives 
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where 

To ;T :T = Gimbal torque commands to the Q ,a gimbal axes 
3xc yyc’ azc ^ x y z 

respectively* 

T :T ;T ^ Desired Control Torques 
cx’ cy C2 

The gimbal command torques are related to the actual generated torque by the 
follovj^ing expression: 
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where 

1/ui^^ = Gimbal torquer tlae constant 

The control torques can have two possible expressions depending 
on the location of the control sensors. For the first configuration 
assume that the sensors are located on body 8 (i.e., the telescope) 
and that the sensor is a rate gyro* For a controller that uses rate, 
position and integral of position the control torque command can be 
written as 
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Substituting equation (2-353) into equation (2-351) yields 
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Substituting equation (2-354) into equation (2-352) gives 
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The second configuration assumes that the sensors are on body 7 • 
Using a similar development , the gimbal control torques can be 
written as: 
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where 


= diag[Kj^,K^y,Kj^] system rate gain 

Kp = diag[Kp^,K^y,Kp^] system position gain 

Kj = diag[Kj^,K^y,K^^] system integral gain 

0)^^ = gyro natural frequency 

? = gyro damping ratio 

s 

2.3 lOG Stability Model - In order to aid in the determination 

of lOG system stability in an efficient and cost effective manner, 

a stability model for the lOG was developed. The stability model 
employed determined the roots of the system characteristic equation 
as a function of parameter variations by determining the eigenvalues 
of the system "A“ matrix. To determine the system ”A" matrix equa- 
tions (2-34) through (2-37) citn be rewritten excluding nonlinear 
terms and disturbances in the following manner: 
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Choosing the system state vector shown in equation (2-362) 
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Eigenvalues of the system matrix as defined by equation 
(2-363) were compared to those obtained from a reduced matrix 
which assumes infinite mass and inertia for the shuttle* This 
reduced *^A** matrix can be essentially obtained by deleting the 
1st and 5th rows and 1st and 5th columns of the matrix de- 
fined in equation (2—363) • The eigenvalues for the reduced "A** 
matrix matched almost exactly with those obtained for the com- 
plete system except for the absence of the roots associated with 
shuttle modes of response. Hence, the reduced **A’’ matrix vja.^ 
used in all of the stability studies performed* 
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Figure 2-3. 
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lOG SYSTEM STABILITY ANALYSIS 


The following sections will describe the lOG system stability 
analyses that were performed during the course of the study. The 
stability studies performed included the following: 

a. System stability as a Function of Telescope Look Angle 
for Nominal System Parameters. 

b. System Stability as a Function of Control Loop Bandwidth 
Variations. 

c. System Stability as a Function of Telescope Mass and In- 
ertia Variations for Nominal Control Loop Bandwidths 

d. System Stability as a Function of Pedestal Mass and 
Inertia Variations. 

e. System Stability as a Function of Variations in Sensor 
and Actuator Characteristics. 

f. System Stability as a Function of Geometrical Variations. 

The nominal system bandwidth that has been used throughout 
this study is 2 Hz defined as the point where the transfer function 
from d .st urbane e torque to command torque is doxjn 3 db. Although 
running higher loop bandwidths would reduce system pointing errors 
due to disturbances originating in the orbiter, it is felt that from 
considerations of system stability in the presence of structural 
flexibility and system noise, using loop bandwidths that are ap-^ 
preciably higher than 2 Hz is not advisable. In addition, the nomi- 
nal suspension parameters used are based on achieving a pointing 
stability of + 1 arc-second in the presence of crew motion dis- 
turbances. The nominal suspension parameters are listed in table 
2-1, A detailed discussion of the stability studies performed Is 
presented in the sections that follow. All of the eigenvalue 
data obtained during this study phase is shown in tho Appendix 
volume which accompanies this report. 

3.1 S ystem Stability as a Function of Telescope Look Angle 
For Nominal System Parameters - For the nominal system as described 
in table 2-1, the telescope was rotated about the y-axis and the 
system eigenvalues were determined at various different look angles. 
The look angle is defined to be zero when the telescope is pointing 
straight up out of the orbiter cargo bay and 90 degrees when the 
telescope is in the cargo bay with its z-axis parallel to the or- 
biter x-axis and is positive directed aft. The telescope look 
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angle was varied between zero degree and 90 degrees for which the 
resulting eigenvalues were all in the left half plan thus verifying 
system stability for the nominal parameters. 

3.2 System Stability as a Function of Control Loop Bandwidth 
Variations - For the nominal system described in table 2-1, the 
pointing control loop bandwidth was varied and the resulting system 
eigenvalues determined as a function of these variations* The 
control loop bandwidth was varied between zero and 10 times the 
nominal loop bandwidth (i.e., 0 and 20 Ez) for telescope look angles 
varying between zero and 90 degrees* The results of these para- 
meterxzations are the following* 

a* All of the system characteristic roots were stable when 
the loop was opened, i.e., the control loop bandwidth was 
zero, 

b. At 1/10 nominal loop bandwidth (i.e*, 0.2 Hz) the system 
was stable for telescope look angles varying between zero and 
40 degrees. For telescope look angles between 60 and 90 de- 
grees, the system was unstable showing an unstable root of 
1*79x10’"^ + j.095 at 60 degrees and 1.48x10“^ j. 1295 at 

90 degrees. These roots describe the interaction between 
translation of the pedestal along the z-axis and telescope 
rotation about the >— axis. 

c. At 1/5 nominal loop bandwidth (i.e., 0.4 Hz) the system 
was stable for telescope look angles between zero and 60 de- 
grees. However, the system was unstable at 90 degrees ex- 
hibiting an eigenvalue of 2.88x10’“^ + j.385. 

d. The system was stable for loop bandwidths varying between 
1/2 and four times (i.e., 1 and 8 Hz) nominal loop bandwidth 
for all telescope look angles between zero and 90 degrees. 

e* Vor loop bandwidths of six times nominal the system was 
unstable for telescope look angles between zero and 10 degrees 
with one set of actuator poles crossing into the right half 
plane. For look angles between 20 and 90 degrees, the system 
was stable* 

f* At ten times nominal loop bandwidth the system was un- 
stable for all telescope look angles with three sets of ac- 
tuator poles crossing to thfe right half plane for telescope 
look angles between zero and 60 degrees and one set of ac- 
tuator poles crossing into the right half plane for a tele- 
scope look angle of 90 degrees. 




3.3 System Stability as a Function of Telescope Mass and 
Inertia Variations - ^or the nominal system described in table 
2-1, system stability was determined as a function of telescope 
mass and inertia variations. The telescope mass and inertia were 
varied concurrently in the same proportion which ranged from one 
hundredth to five times the nominal values shown in table 2-1, 

During these variations the pointing control loop bandwidth was 
continually adjusted to maintain a control loop bandwidth of 2 Hz. 
All other parameters were maintained at their nominal values. The 
results indicate that the system was stable for all values of tele- 
scope mass and inertia considered for telescope look angle ranging 
between zero and 60 degrees, 

3.4 System Stability as a Function of Pedestal Mass and In- 
ertia Variations - For the nominal system described in table 2-1, 
system stability was determined as a function of lOG pedestal and 
inertia variations. The pedestal mass and inertia were varied 
concurrently in the same proportion which ranged between one tenth 
to 10 times the nominal values indicated in table 2-1. These vari- 
ations were conducted for zero and 60 degree telescope look angles. 
The results of these studies indicated that the system was stable 
for all variations made. 

3.5 System Stability as a Function of Sensor and Actuator 
Variations - For the nominal system described in table 2-1, system 
stability was determined as a function variation in sensor and ac- 
tuator characteristics. The gimbal torquer (modeled as a first 
order lag) time constants were varied from 20 to two times their 
nominal values (i.e. , from .1591 to 1.591xl0"2 sec). The rate 
gyros (which were modeled as second order transfer functions) 
natural frequencies were varied from one one-twentieth to one half 
times their nominal values (i.e., from 8.17 to 81.7 rad/sec) and 
damping ratio was varied from one hundredth to 0.6 times their 
nominal values (i.e,, 0.005 to 0.3). The results of these studies 
were the following: 

a. The system was stable for the range of gimbal torquer 
variations investigated for telescope look angles between 
zero and sixty degrees, 

b. System was unstable when the gyro damping ratio was 0.01 
or below and was stable when the sensor damping was 0.05 (i.e., 
one-tenth nominal) or greater for telescope look angles be- 
tween zero and jixty degrees. 


j 
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c. System was unstable when the gyro natural frequency was 
16.34 rad/sec or below and was stable for gyro natural fre- 
quencies of 32.68 rad/sec or greater for telescope look- 
angles varying between zero and 60 degrees. 

These results indicate that system stability is not a sensi- 
tive function of variations in sensor and actuator characteristics. 

3* 6 System Stability as a Function of Variations in Suspen- 
sion Characteristics - For the nominal system described in table 
2-1 , system stability was determined as a function of variations 
in suspension damping and linear stiffness. The suspension damp- 
ing and stiffness were first varied together from one-tenth to 
20 times the nominal values shown in table 2-1. Then with the 
suspension stiffness set at nominal values, suspension damping 
was varied from zero to one-half nominal. The results of these 
studies are outlined below; 

a. The system was stable when suspension stiffness and damp- 
ing were varied together over the range indicated above for 
telescope look angles between zero and sixty degrees. 

b. The system was unstable for all telescope look angles 
when suspension damping was set to zero. 

c. With suspension damping set at one-tiundredLh the nominal 
values shown in table 2-1 keupxng the stiffness at their 
nominal values, the system was stable for telescope look 
angles between zero and 60 degrees. For a ninety degree 
telescope look angle, the system was unstable. The above was 
also true when suspension damping was set at one-fiftieth 
nominal. 

d. The system was stable for telescope look angles between 
zero and 90 degrees when suspension damping was one-tenth 
nominal or greater. 

These results indicate that systeta stability is not a sensi- 
tive function of variations in suspension characteristics. 

3.7 System Stability as a Function of Geometrical Variations 
For the nominal system described in table 2-1, system stability 
was determined as a function of geometrical variations. These 
variations included the following: 

a. Telescope CM was moved in the x,y plane, (i.e., along 

the telescope transverse axes) between the nominal value of 
zero to two meters along the x and y axes. 
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b. The distance of the telescope CM with respect to the 
gimbal hinge point, r 2 > was varied from the nominal value of 
1.689 meters to 168.9 meters. 

c. The location of the gimbal hinge point with respect to the 
pedestal CM, r^, was varied along the x and y pedestal axes 
between the nominal value of zero to 0.375 meters, 

d. The variations described in items a and c above, were per- 
formed concurrently. 

The results of these studies are described below: 

a. When the telescope CM was moved along the transverse axes 
of the telescope (i,e., x and y axes) the system was stable 
for offset values along the x and y axes of as much as two 
meters for telescope look angles varying between zero and 
60 degrees* 

b* When the distance of the telescope CM with respect to the 
gimbal hinge point, r2» was 84.55 meters or less, the system 
was stable for zero telescope look angle. However, for tele- 
scope look angles of 60 degrees, the distance of the telescope 
CM with respect to the gimbal hinge point had to be kept to 
33.78 meters or below for the system to be stable. 

c. The system was stable for variations in the location of 
the gimbal hinge point with respect to the CM of the pedestal 
r]^, along the x and y axes of as much as 0.375 meters which 
corresponds to the distance that the gimbal hinge point is 
above the pedestal CM (i.e., distance along the z-axis) for 
telescope look angles between zero and 60 degrees* 

d. The system was unstable when the concurrent values of 
r^ and r 2 were 
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for a telescope look angle of zero degrees. When the telescope 
look angle was 60 degrees, the system was stable for the above 
conditions • 

These results indicate that system stability is not very sensi- 
tive to realistic geometric variations that will reasonably be 
encountered. 

A summary of the results obtained in the lOG stability studies 
performed and described above is given in table 3-"i. 
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Table 3-1. Summary of lOG Stability Studies 


study Performed 

Parameters Varied 

Range of Variation 

Results 

System stability as a function 

Suspension linear (D) 

0<D< 0 . 5 nominal 

• For zero damping coefficient 

of suspension damping for var- 

and angular (d) damp- 


system was unstable for all 

ious telescope look angles 

ing coefficient 

0£d<^ 0.5 nominal 

telescope look angles. 

• For 1/100 nominal damping system 

2 Hz control loop bandwidth 

Telescope look angle (6) 

_ 

— 

O<0<_ 90“ 

was stable for look angles of 
60® or less. 

• Unstable for telescope look 
angle of 90® . 

• System stable for all telescope 
look angles for 1/10 nominal 
suspension damping 

System stability as a function 
of suspension characteristics 

Suspension linear (K) 
and angular (k) stiff- 

0.1 nom ^K<20 nom 

•All conditions examined were stable. 

keeping damping ratio constant.j 

ness. 

0.1 nom £k^20 nom 


Telescope look angle 60® 

Suspension linear (D) 
and angular (d) damp- 

0.1 nom £D^20 nom 


|2 Hz control loop bandwidth 

ing coefficient 

0.1 nom nom 


System stability as a function 
of pedestal mass and inertia 

Pedestal mass (m 2 ) 

0.1 nom 

•All conditions examined were 
stable. 

variations 

Telescope look angle 60® 

2 Hz control loop bandwidth 

Pedestal inertia (J 2 ) 

0.1 nom 
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Table 3-1. Summary of lOG Stability Studies (Continued) 


Study Performed 


of sensor damping ratio for 
various telescope look angles. 


2 Hz control loop bandwidth 


System stability as a functi 
of sensor natural frequency 
for various telescope look 
angles . 

Sensor damping ratio 0.5 
2 Hz control loop bandwidth 


of actuator band width for 
various telescope look angles. 

2 Hz control loop bandwidth 


Parameters Varied 

Range of Variation 

Results 

Sensor damping ratio 
C?s) 

Telescope look angle 

O) 

0*01 nom <£ < nom 
O<0£ 60“ 

• System unstable for sensor damp- 
ing ratio less than .05 which is 
1/10 the nominal damping ratio. 
This was true for all look angles 
; considered. 

Sensor natural fre- 
quency (W^g) 

0.05 nom nom 

— ns— 

• System unstable for sensor natural 
frequency less than 1/5 nominal 
(i.e. 32.68 rad/sec). This was 
true for all telescope look angles 

Telescope look angle (6) 

o<e< 60“ 

considered. 

Actuator bandwidth 
Telescope look angle (6) 

0.05 nom tiom 

O 50 < 60“ 

•All conditions examined were 
stable. 













Studies (Continued) 


Range of Variation 


Results 


-.375<rj^j^£ .375 


-.375<r, < .375 
“ ly- 


r, =.375 (nom) 
Iz 


-2.0<r, <0 
— 2x— 


0<,t2y~^ • 


1.689<r22£l68.9 


0<6< 60' 


•System stable for separate vari- 


ations in for all conditions 


considered , 

•There were only two conditions 
examined for this parameterization 


= -.375i +.375i +.375i_;e=0%60^ 

1 X V 2 


,375i3^-.375iy+.375i^; 6=0' 


* System was unstable for 
= .375i2(nora) 


r^<85i ;0=O 
2— z’ 


r,<33.78i ;0=6O‘ 
2— z’ 


• System was stable for nominal r^^ 
and nominal X2z 1*689) 

for the total range of parameter- 
ization of r_ , r„ and 0. 


2x’ ^2y 

•System was unstable for the follow- 
ing condition 


r, = .375i -.375i +.375i 
1 . X ^ y . z 


= 1.5i H-1.5i,>l,689i 

2 X y z 


e = 0" 

When the telescope angle was equal 
to 60® the system was stable for 
the above geometric conditions. 




System stability as a function Telescope mass (m 2 ) 
of telescope mass and inertia 
variations for various 

telescope look angles Telescope inertia (J 2 ) 


Telescope look angle 6 


Studies (Concluded) 



Range of Variation 

Results 

0<f < 20 Hz 
“ n— 

•System was unstable for control 
loop bandwidths of .2Hz (1/10 

0<6< 90“ 

nominal) or below when telescope 


look angles were 60 and 90 degrees. 

•System was unstable for control 
loop bandwidths ,4Hz (1/5 nominal) 
or below when the telescope look 
angle was 90 degrees . 

•System was stable for all telescope 
look angles for lHz<f <4Hz and un- 
stable for f pg6Hz . 

. 01 nom ^m2<5 nom 

.01 nom ^J2;<5 nom 
O<0£ 60“ 

•All conditions eKamined were 
stable. 
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lOG LINEAR SYSTEM PERFORMANCE 


The sections that follow describe the pointing performance 
of the lOG for nominal conditions and as a function of various 
system parameterizations which include the following; 

a* Location of the lOG with respect to the orb iter* 

b* Variation of suspension damping characteristics. 

c. Variation of pedestal mass and inertia. 

d. Variation of telescope mass and Inertia. 

In the performance of the system parameterization studies 
outlined below, a nominal set of system parameters were employed 
which are given in table 2-1. These parameters result in a 2Hz 
pointing control loop bandwidth defined as the point where the 
magnitude of the ratio of control torque to disturbance torque 
is down 3 db. Any deviations from these parameters will be ex- 
pressedly noted in the discussion that follows. The crew motion 
disturbance force profile that was employed throughout the point- 
ing performance studies is shown in figure 4'“1. 

lOG nonlinear system performance will be described in section 

5.0. 


Earth tracking performance of the lOG for the nominal system 
given in table 2-1 will also be described in the sections that 
follow. 

4.1 IQG Pointing Performance 

4.1.1 IQG Pointing Performance as a Function of Location With- 
in the Orbiter Cargo Bay - The nominal parameters given in table 
2-1 position the lOG as far forward in the orbiter cargo bay as 
possible. In order to position the lOG as far back in the cargo 
bay as is feasible, the baseline parameters shown in table 2-1 are 
modified in the following manner, 

- 4.18 1 - 1*579 1 m 

120 X 2 

System responses under the influence of crew motion dis- 
turbances, applied along all three vehicle axes simultaneously, 
were taken for various telescope look angles for both positions 
of the lOG within the cargo bay. The peak pointing error incurred 
for the X and y telescope axes as a function of telescope look 
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angle for both lOG positions within the cargo bay are shown in 
figures 4-2 and 4-3, Examination of these figures indicate that 
pointing errors incurred when the 106 is mounted forx-jard in the 
cargo bay far exceed those that result when the lOG is mounted 
in the rear of the cargo bay. The reason for this behavior is 
the following. When the lOG is mounted in a forward position the 
orbiter rotations and translations that occur due to the crew 
motion disturbance aJd to give the net linear motion of the lOG 
hinge point which is the primary disturbance input to the tele- 
scope. However, when the lOG is mounted in the rear of the cargo 
bay the orbiter translatory and rotary motions subtract to give the 
net translation of the of the gimbal hinge point thus resulting in 
a greatly reduced disturbance input to the telescope which in turn 
results in small pointing errors. It is for this reason that the 
baseline system parameters were defined so that the lOG Is in a 
forward position in all of the ensuing response studies. 

4.X.2 IQG Pointing Performance as a Function of Variation in 
Suspension Damping - System time response studies were performed 
for various telescope look angles as a function of suspension damp- 
ing characteristics. The suspension damping was varied over a wide 
range of values and the effect on peak pointing error for a reduc- 
tion in suspension damping of an order of magnitude is summarized 
in figure 4-4. The parameter values that result in a factor of 
ten reduction in suspension damping are given below. 

D = D = D = 5.993 — 

X y z m 


d = d = 1,496 n-m-sec 
X y 

d = 2.992 n-m-sec 
z 

Examination of figure 4-4 indicates that suspension damping has 
very little effect on system pointing performance when varied (i.e., 
reduced) by an order of magnitude. However, suspension damping 
should not be reduced by much more than an order of magnitude if 
adequate system stability margins are to be maintained. In addi- 
tion, the linear mode of system oscillation is not damped by the 
control loop and only damps out by virtue of the damping in the 
suspension system. The nominal suspension damping coefficients 
are chosen such that a damping ratio of a tenth is achieved for 
the linear mode when considering all of the mass above the sus- 
pension point. 
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Figures 4-5 thru 4-25 show the lOG system response for nomi- 
nal parameters in the presence of tht* crew motion disturbance 
shown in figure 4-1 applied simultav.eously along the x,y, and z 
axes of the shuttle orbiter, for a telescope look angle of zero 
degrees* The legend that is employed in these plots and all sub- 
sequent plots sho™ in this report is given in table 4-1. 

The symbols employed in table 4-1 are the same as those used 
in the mathematical model derivation outlined in section 2*1 and 
are defined in that section. It should be noted that v;hen 

and are equal to zero, which they are in the case of stellar 
sV it 

pointing, and are identically equal to cjt)^ and 0^ the angular 

rates and rotations of body 3. Figures 4-26 thru 4-46 show the lOG 
system response for nominal parameters with a crew motion disturb- 
ance applied simultaneously along all three orbiter axes for a tele- 
scope look angle of 0.6981 rad (40 degrees) while figures 4-47 thru 
4-67 show lOG response for the same set of conditions but with a 
telescope look angle of 1.047 rad (60 degrees). 


4.1.3 lOG Pointing Performance as a Function of Pedestal 
Mass and Inertia Variations - Time response studies were per- 
formed for various telescope look angles as a function of pedes- 
tal mass and inertia variation. The pedestal mass and Inertia 
was varied by an order of magnitude about nominal (i.e. , in both 
directions) with both the mass and inertia varying concurrently 
by the same ratio from nominal. The parameters that were modi- 
fied from the nominal set shown in table 2—1 in order to accomp- 
lish the desired pedestal mars and inertia variation are given 
below: 


- 19.5 Kg 




One tenth nominal 
pedestal mass and inertia 


m2 ^ 1950 


J 


2x 



Kg 

= = 500 Kg-m^ 

Zz 


Ten times nominal 
pedestal mass and inertia 


The results obtained from these time response studies are 
summarized in figure 4-68. Examination of this figure indicates 
that variations in pedestal mass and inertia by an order of magni* 
tude from nominal docs not have a marked effect on pointing per- 
formance. It can be seen that when the pedestal mass and inertia 
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is reduced by au order of magnitude the x-axis pointing performance 
improves while the y-axis performance is unchanged* When the ped- 
estal mass and inertia are increased by an order of magnitude the 
K-axis pointing errors increase while the y-axis errors decrease* 
This is to be expected since the suspension does not act primarily 
as a classical passive isolator in reducing errors in the Xj y 
plane. It accomplishes that function by allowing the hinge point 
to translate ”via" pedestal rotation about an appropriate axis 
in the x,y plane whose center of rotation is between the pedestal 
CM and the plane defined by the attach point of the individual 
isolators to the pedestal base plate. The rotation of the pedestal 
takes place due to the geometrical properties of the lOG and the 
influence of the telescope control torques which react back on 
the pedestal* Hence, as the pedestal mass and inertia decrease 
the speed with which the pedestal rotation in the x,y plane occurs 
with its accompanying hinge point translation increases thus re- 
ducing the disturbance input to the telescope, which results in 
improved pointing accuracy* When the pedestal mass and inertia 
are increased, the rotation of the pedestal in the x,y plane 
occurs at a slower rate which results in larger net hinge point 
accelerations thus increasing the disturbance input to the tele- 
scope which results in larger system pointing errors* For y-axii? 
disturbances, the lOG suspension acts primarily as a classical 
passive isolator* Hence as the pedestal mass and inertia is 
increased the system pointing performance improves since the 
isolation natural frequency is lowered. However, it is seen that 
there was virtually no effect on the y-axis pointing error 
when the pedestal mass and inertia was reduced by a factor of ten. 
The reason for this is that the mass that is used to compute the 
linear suspension natural frequency in the **y" direction is a 
composite mass that is a function of the pedestal and telescope 
masses and the pointing control loop bandwidth. For the nominal 
pedestal mass and inertia the composite mass which determines 
suspension natural frequency is primarily determined by the tele- 
scope mass and the pointing control loop bandwidth. Hence re- 
ducing pedestal mass and inertia has virtually no effect on 
the y axis pointing error* This effect will be further demon- 
strated by the results obtained when the suspension linear stiff- 
nesses in the x and y directions are increased while the rotational 
stiffness about the x and y axes are held constant. 

4*1*4 IQG Pointing Performance as a Function of Telescope Mass 
and Inertia Variations - Pointing performance studies were performed 
using a small telescope defined in Ddrnier technical report IPS— TN- 
025-041. The mass properties of this telescope given in the above 
report are 

m = 290 Kg 
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= Jy = 50 Kg-m^ 

= 35 Kg-m^ 

• -A 

Distance from the Base of Telescope to C.M. ■‘1.5 ^ 

lifhen these parameters are combined with the lOG stable mount- 
ing gimbal the following are the parameters tliat have been modified 
from the nominal set shown in Table 2-1. 

m^ - 583 Kg 

J3^= 309.5 Kg-ra^ 

J^y = 334.5 Kg-m^ 

J_ = 110 Kg-ra^' 

3z 

■ A.' • 

= 0.9346 i m 

A z 

= diag [3.032x10^; 3.277x10^; 1.078x10^] 

' Kj, = diag [1.076x10^; 1,163x10^; 3.824x10^] 

Kj. = diag [1.698x10^; 1.835x10^; 6.035x10^] 

where a 2 Hz pointing control loop bandxiridth was maintained 
for the small telescope. 

The results of the pointing performance studies are summar- 
ized in figure 4-69, Examination of figure 4-69 Indicates that the 
pointing error has increased by approximately a factor of five 
over: those for the baseline telescope for both the x and y axes. 
This is not surprising since the loop bandwidth was maintained 
at 2 Hz for the small telescope and hence one might expect that 
the pointing error for the small telescope could be approximately 
derived from those Incurred for the baseline telescope using tlie 
following relationship. 


n-m-sec 

n-m/rad 

n-m 

sec 


®cils 


"2's 


2'nom 


where 




3i ' nom 


'31 's 




ei'npm 


(4-1) 


~ Distance from hinge point to telescope CM for small 
telescope 

r„| = Distance from hinge point to telescope CM for base— 

2 line telescope 


J„. I = Moment of Inertia of the nominal telescope plus the 
X nom JOG stable mounting girahal about the 1^^ axis i = x,y 

= Moment of Inertia of the small telescope plus lOG 
stable mounting gimbal about the i*^^ axis 1 = x,y 


0 . I = Pointing error incurred by the nominal telescope 
X nom the i*-^ axis i = x,y 


0p. I = Pointing error incurred by the small. telescope about 
^ ® the i*^l^ axis i " x,y 

Substituting the appropriate parameter values in equation (4-1) 
and evaluating results in 


0 

~ 5.190^ 

ex's 

Ex 'nom 

0 1 

= 4.380^ I 

ey'B 

ey'nom 


(4-2) 


It is seen that equation (4-2) gives approximately a factor of 
five increase in pointing error for the small telescope over that 
incurred by the large telescope which is in good agreement with the ; 
simulation results. 

In order to achieve performance similar to that achieved for 
the baseline telescope Ci.e, , between 0. S and 0. 6 arc-seebtid) the 
suspension stiffness had to be reduced by a factor of five. These 
results indicate the strong possibility that one set of lOG sus- 
pension parameters will not yield sub arc-second ppinting perform- 
ance over the spectrum of payloads that are being considered for 
lOG stabilization . It is probable that two or perhaps three dif- 
ferent suspension system charncteristics xvill be required to yield 
sub arc-second pointing and tracking performance over tlie payload 
set being considered. 


4 
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4.2 lOG Earth Tracking Performance -- The earth point track- 
ing performance of the baseline system defined in table 2-1 was 
determined for a 438 KM orbit. Both angular rate and position 
commands to the telescope pointing control system were used simul- 
taneously to accomplish earth tracking* These profiles are given 
in equations (4-3) and (4-4) * 


(4-3) 


(4-4) 

The profiles shox-m in equations (4-3) and (4-4) are actually 
approximations to those derived in section 8> however, they are 
"good** approximations as indicated in that section and hence yield 
representative system performance when they are employed. 

The results of the earth tracking studies for the nominal 
system indicat-e that the system tracking error was 0.86 arc- 
second peak which is x^ithin the desired tracking performance of 
41 arc-second \j±th an accompanying peak control torque of 5.33 
n^m which is within the 20 n-m gimbal torquer capability. However, 
the pedestal rotation about the y axis 0„ x^as 0.249 rad (14.27, 

deg) and the pedestal translations were 9.92x10 m (3.906 in) 

■-3 ■ 

and 4.39x10 (.1728 in) along the x and y axes respectively. 

These values resulted in isolator elongations of 2.60 in along 
their longitudinal axes. These values of pedestal rotations, 
translations and isolator elongations are x-jell beyond the lOG 
suspension design limits and hence the nominal suspension sys- 
tem cannot accommodate the nominal payload for earth tracking from 
a 438 KM orbit. It should be noted that the situation becomes 
worse as the orbital altitude is decreased (see section 9) and 
the ensuing pedestal rotation translations and isolator elongations 
will increase further^ It is, therefore, evident that if the lOG 
is to accommodate payloads in the 200 Kg class the suspension 
parameters must be modified, consis taut with painting performance 
requirements , to reduce pedestal rotations and translations. A 
possible approach to this problem i 
pension linear stiffness along the 
the rotational stiffness about the 
stiffness along the z axis at the i 
The rationale behind this approach 
act as a classical isolator along 
ing the linear stiffness along thei 
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grould be to increase the sus-^ 

X and y axes and to increase 
X and y axes keeping the linear 
lominal value shovm in table 2-1. 
is that the isolator does not 
the X and y axes, hence, increas- 
se axes will have little effect 


0^(5) = 1.136 tanh[-5||-A.] 

0) (t) = -1 . 748xl0~^sech^ 
c. 65 


on the telescope pointing error incurred about these axes due to 
disturbance as originating in the orbiter (e,g*., crew motion) . 
Increasing the rotational stiffness about the x and y axes will 
increase the pointing error incurred about tbe x axis due to or*- 
biter disturances, however, the pointing errors presently obtained 
with the nominal system about the x axis Is small enough (i.e«, 

0.1 arc-^second) to allow an increase in this error while still 
meetitig the overall pointing stability requirement of 4?- arc- 
second. A more detailed discussion of this subject and the re- 
sults of the approach outlined above will be given in section 
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Figure A-13* Pedestal z-Axte Rotation 
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Figure 4-44* x-Axis Control Torque 


TELESCOPE LOOK ANCLE * 40* 

CREW MOTION DISTOHDAHCE ALONG THREE AXES 







RIJN ?l '4 


Figure 4-47. 
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Figure 4-68o Peak Pointing Error as a Function of Pedestal Mass and Inertia Variations 
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lOG LIMIT CYCLE ANALYSIS MD NONLINEAR FERFOIttlANCE 


The sections that follow describe a limit cycle analysis 
using describing function techniques for the lOG, under the 
influence of gimbal wire torques. Results obtained from this 
analysis are presented and compared to those obtained from the 
three axis lOG time simulation which gives a better understand- 
ing of the limit cycle phenomena observed and provides a cross- 
check to the time simulation. 

The three body Shuttle* lOG pedestal, lOG inner gim- 

bal/ telescope combination) three axis lOG equations of motion 
described in section 2,1 were reduced to a tt-^o body system con- 
sisting of the lOG pedestal and lOG gimbal/ telescope combination. 
The analysis technique presented is three axis handling the three 
gimbal wire torque nonlinearities simultaneously for arbitrary 
telescope look angles* In order to facilitate numerical compu- 
tation these equations were subsequently simplified by only con- 
sidering planar motion and gimbal wire torques about the y axis 
of the lOG. The planar motions considered had the following 
degrees-of-freedom: 1) Rotation about the lOG pedestal y-axia; 

2) rotation about the y axis of the lOG inner gimbal/telescope 
combination; 3) pedestal translation along the x axis; 4) pedestal 
translation along the z axis. lOG system limit cycle behavior 
was quantitatively investigated for telescope look angles of zero 
and 40 degrees and subsequently compared to results obtained from 
the time simulation. In addition, time response plots for various 
wire torque characteristics for a telescope look angle of 40 
degrees are presented* 

5.1 Describing Function Analysis of IQG System 

5*1.1 IQG Wire Torque Characteristics - The wire torque 
characteristic employed in this analysis shovm in figure 2—2 is 
repeated here for convenience. 

As described in the December 5, 1975 bimonthly report of 
NASA Contract Number NAS8-31570 titled, "Research Study on IPS 
Digital Controller Design," the i/ire torque characteristic shoira 
in figure 5-1 can be reduced as shown in figure 5-2. 

The form of the wire torque characteristic shoxra in figure 
5-2 is easier to handle analytically and will be used throughout 
this analysis. 




5.1.2 lOG Equations of Motion - The equations of motion for 
the lOG, considering the system as three bodies consisting of 
the shuttle, 106 pedestal, and lOG inner gimbal/telescope assem- 
bly are derived in section 2.1. The fina], equations of motion 
given in section 2.1 are repeated here for convenience. 
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where 


^WT^2 ~ torques written in body 2 coordinates 


^WT ' 3 ~ torques written in body 3 coordinates 


,T„ - „T„ 


-1 


= Trans Cormation between bodies two and three 


I 


The "12** subscript used in the derivation sho™ in section 2*1 
is dropped here, for convenience. 

Since the shuttle mass and inertia far escceeds that of the 
lOG and the telescope, the above equations can be simplified by 
taking the limit as the shuttle mass and inertia {i,e,, and 
approach infinity. This results in the following set of equations 

’^2e+®22-^2e'^l-2^3-^3e-™WTl2 ’ '‘^0-'”’^20‘“2«’®E20-®2> 
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^2e**'2^3‘ ^3e ^*'^®*^20*“2~^^‘®E20*®2'*“3^1*‘^2 




(5-7) 


vjhere the nonlinear terms shown in equations (5-2) through (5“4) 
eliminated o 


In order to make use of the wire torque characteristics shown 
in figure 5-2, the equation of motion for the lOG must be written in 

terms of relative accelerations, rates, and angle between the pedestal 
and the telescope (i,e«, between bodies 2 and 3), This is accomp- 
lished by making the following substitution 

“3 ' 3'^2'“2-^ 


“3 ' 3^2 -“2+r 


«3 ' sW^ 
®3 ' 3’^2’«h^ 


r\ 


Substituting equation (5-8) into equations (5-5) through (5-7) gives 

T2e+“22’^2a-^l-2''3-^3e'™Hll2 " ^•'J-“3'l‘2’^3''2‘3’^2^'“2 

+(d-Alij22o*®'®U20''2S**Sl*3'^2^*“2'^^*^”^*^20’^‘^20"2^3**^P*3^2^*®2 

-t2l3-Ki-3T2)-9^tM3?,-2T3-72l-Y-2T3-!Kj-Y+Kp-T^Kj-7b 

. e+4Rp20 ‘ ^ 1 
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+ (m 2 +m 3 ) e+4D • £+4K* e 

• «■ 

The relative angle, rate, and acceleration (i.e, Y, Y» Y) 
are written in telescope (body "3”) coordinates. However, these 
quantities should be written in terms of the relative Euler angles 
between the lOG pedestal and telescope (i.e,, bodies ”2" and ”3”) 
since this will be required if single input describing functions 
are to be employed in the description of the gimbal wire torques. 


The rotation sequence between the lOG pedestal and the tele- 
scope is 2 , X, y. Defining these rotations as 0 , respect- 

ively, the relative rate Y can be written as; 
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As a consequence of the nonlinearity breakdown shown in sec' 
tion. 5.2.1 the gimbal wire torques T- , , can be written as 
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where 


N = the describing function of the switch shown in figure 2 
= the wire torque slope (linear) shown in figure 2 

and 


(5-19) 

( 5 - 20 ) 


N = dlag [Ng^, 

= diag [Ky^g^s 

Substituting equations (5-13), (5-15), (5-17) and (5-18) into equa- 
tions (5-9) through (5-11) the following results 
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V2''3'''3a=-<Vl^32’'3‘^2’3'f2’-“2-“'®E20'“2-«'®E20-®2-”32^3‘’^2-^‘^ 
+(m 2 *ha 2 )^ +4D*£+4K*e (5-23) 

It should be noted that in the derivation of equations (5-21) 
through (5-23) the "A” and "B” matrices were considered to be 
constant. In fact " 3 ^ 2 ” ^i^d " 2 ^ 3 " also considered constant 

if the system is to be linear where a describing function analysis 
on the gimbnl torque nonlinearity can be conducted. These appi'oxl- 
tnations are valid since the limit cycle amplitude Incurred Is 
generally **sm.ill" and the transformations can be considered con- 
stant. 
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The system equations derived above (excluding disturbances) can 
be written in the following form 
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The system matrix given In equation (5-27) can be used to 
generate three axis root loci as a function of variations in the 
describing function amplitude K. The condition for a limit cycle 
to exist will be for those values of N for which there are poles 
on the "juj” axis. The stability of the limit cycle can be easily 
established by noting whether the root locus crosses from the left 
half to the right lialf plane or vice-versa as a function of de- 
scribing function amplitude variation in the vicinity of a parti- 
cular pole on the axis. If the locus goes from the left 

hand plane to right half as the describing function is Increased 
the limit cycle is stable if the opposite occurs then the limit 
cycle is unstable. It should be noted, however, that in the above 


(5-28) 


5-9 


derivation it was implicitly assumed that a single sinusoidal 
input is sufficient to characterize the input to each of the 
glmbalwire torque nonlinearities r This assumption has proved 
to be adequate when analyzing 10'; planar motion as described in 
the sections that follow. The adequacy of the assumption when 
all three girobal nonlinearities are present with the telescope 
at an arbitrary angle with respect to the lOG pedestal which 
cross couples the lOG axes still remains to be demonstrated. 

This demonstration will come when the analytical results obtained 
from the root locus approach described above are compared with 
simulation results obtained from the time simulation. 


Examination of equations (5~21) through (5-23) indicates that the 
limit cycle characteristics do not depend on since it is 

generally much smaller than K^. 

5.1,3 Planar Limit Cycle Analysis - Equations (5-21) through 
(5-23) which describe three dimensional motion of the 100 can be 
reduced to planar motion which consider the following degrees of 
freedom t *)2 , and These simplified equations are 

listed below. 
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I 




ii ' 

i f 


'i; 


r" 

\ ; 
; 




Taking the laplace transform of equations (5-29) through (5-32) gives the 
folloT-jlng: 

W®22’''2jy-^l-2’^3‘^3Jy ' <<''2y-^3’^L'^3’^lzV‘’='<'y<.^®^'^<^y’““*®^20z-'4ly^° 

K K 

+V“*®?20z-V>- -f ^®2y-^"3V22“=V'- “f ^ 

(5-33) 

K 

’^3zy+'^3'^2>'''3aly = ^<-’3y-^34z'^3V22“='l'y<>’®^'^‘'Ky®^Py’^ 

(5-34) 

W2’^3"^3elz ‘ {("3’^lz‘*”3’^2z“='''yo’=^‘”\“2202=+«x'*E202’®2y 
'^3^2z‘^°®’^yo^^'^y^^ (^2*^3^ (5-35) 

^2ezV3‘^2elz " -"3"2z®‘“'''yo^^®2y-”3'2z=^'‘V®S*‘'"2'^3^®^'*^V+“z>^2 

The system characteristic equation is given by 

-N 


N 


A = det 


(5-37) 


r- 


i: 
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Expanding equation (5-37) yields 


A = a? ag 


^9 -^^9 ^4 


which can be fui:th^?r expanded to 

A « 02 [a^ag+ag (Ogag-a^ag) l+a^ag (aj^ag-agag)+aga^(ag*-ag)+agag (aga^-aj^ag)-a ^^a^a 

(5-39 


K 

' '•’2y-^3’^L^3’^lz-^2z“=V>='-"V'’”x420-V>='^<V“x'4o-Sy)- -f 

“2 “ Vlx''2x“=^o'^‘%'™i|y^®‘<V‘Wy> " 

=3 ° Vu^^'”®E20zV-^®E202‘'x 
% “ «-’3y-^3'L-^3"lz’^2. “=V>'Hy"+V o 


4 = W3y'^3’^L>='+< W^^'^^py'^y’-' 


(5-4C 


a, = in-,r„ cos'/' S'^ 
6 3 2z ^yo 


a- = -m„r sini(» 

7 3 xy ’^yo 


ag - (m3r^^-Hn3r2^cos'/' )S +^\Re20z^^'‘^x^E20z 


a„ - (m„-hn^)S +4D S-I-4K 








For the system presently under consideration the following parameter 
values apply: 


IQn 




« 195 Kg 

= 2.683x 103 Kg 
2 


‘'2y 

■^3y 

^Iz 

^^2z 

^202 

d 


3 „2 


= 50 Kg-ra' 

« 2.648x10-^ Kg-m 
= .375 m 
•» 1,689 m 
= -.375 m 

n-mrsec 


D 


= 14,96 
= 31.25 
= 59.93 


rad 

n-m 

rad 

n-sec 


m 


K = 125 - 

X m 


% 




= 1.972x10 


4 n-m-sec 


rad 


“ 7.00x10 


*=* 1.105x10 


4 n-m 
rad 

5 n-m 
sec 


5. 1.3.1 Limit Cycles for Zero Telescope Look Angles - Assuming 
*Sro|^y ” '^yo ” ^ following results upon substitution 

of the system parameter values into equation (5-40), 
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I 

i; • 


f' ‘ 


2. 12664x10^8^-1. 96713x10^8-6. 98984x10^- 


a” = 1. 69935x10^8^- (1.972xl0^+Nj^y^S- (7. 0025x10^)- 


.3„2 


= 1.00613xl0‘'S -89.895S-187.5 


(5-41) 


= 1. 20012x10'^ S^+(l. 972x10^) SH~(7 . 00x10^ )■)- ^ ‘ ■ 


V 


a” = l.O3O19xl0^S^+(1.972xl0'^+N, ) S+(7. 0025x10^)+ 

D yy b 


a, = 4.53159x10^8^ 

D 


- 0 


3 2 


ag *= 5.53771x10^8 -89.8958-187.5 


^3„2, 


ag « 2.878x10 S +239.728+500 


Substituting equation (5-41) into equation (5-39) and evaluating the 
expression yields: 


A = 6 .5039xl0®S^+l. 5783xl0^°8^+7 . 283xl0^‘^S^+2 . 2356x10^^8^3-57909x10^^8^ 


+2.87111x10^^8^+2.6x10^8+1. 7x10 Vn, [9.99371x10^8^4.8795x10^8^+9.437x10^8^ 

'i'y 


+1.4x10^8+1.56x10^] = 0 


(5-42) 




Rearranging the terms in equation (5-42) gives 


8^ [9 . 99371xl0^s'^+4. 8795x10^8^+9 . 437x10^8^+1 . 4x10^8+1 . 56x10^ 1 


(6. 5039x10^8^+1. 5783x10^^8^+7. 283x10^^8^+2. 2356x10^^8 V3. 57909x10^^8^ 


N,. 


+2 . 87111x10^^8^+2 . 6x10^8+1 . 7x10® ) 


>7 

(5-43) 


Substituting S=j0J into equation (5-43) results in 

CJ^{ (9. 99371xl0^a>^-9.437xl0^m^+l. 56xl0^)+1 (-4. 8 795x10^(0^+1. 4xl0^‘g)1 > 


[-1 . 5783xl0^‘^w^+2 . 235 6x10^ ^ uj‘^-2 . 87111xlO^\s^+l . 7x10®] 


N, 


+j (-6 . 5039xl0^u7+7 . 283xl0^°ui^-3 . 57 909x10^^ w^+2. 6xlO®oi] 


l|»y 

(5-44) 


1 
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In order for equation (5-44) to be valid the imaginary part 
must be equal to zero. Rationalizing equation (5—44) results in 
the following necessary condition, 

in 9ft 9A 9 

w "1.01075x10 w +4.88179x10 w "3.07797x10 (*) +2.73930) "2,57854x10 = 0 

Solving equation (5-45) for w gives 

W = +.03271 

U) = +.08914 

to — +.8564 

to = +2.0756 

to = +9.7992 

In order to determine whether all of the solutions for tt) 
given in equation (5-46) are possible limit cycle points, the various 
solutions will have to be inserted into the real part of equation 
(5-44), given in equation (5-47), the results evaluated and shown to 
be positive 

Real Part = 0)^{ (9. 99371x10^0)^-9. 437xlO®to^+l. 56x10^) 

(“1 . 5783xl0^°U)^+2 . 2356xl0^^0)'‘^-2 . 87111x10^^0)^+1 . 7x10^) 
+(-4. 8795x10^0)^+1 . 4x10 V) 

(-6 . 5039x10^0)^+7 . 283xlO^°0)^-3 . 57909x10^^0)^+2 . 6x10^0)) } 

(-1. 5783xl0^°o)^+2 . 2356x10^ V-2 .87111x10^^0)^+1 . 7x10^) ^+(-6 . 5039x10 V 
+7 . 283x10^ V-3. 5 7909x10^^0)^+2 . 6x10% 

When the above procedure is followed the solution of 01=9,7992 rad/sec 
is not a possible limit cycle point while all the others are. . 

The amplitude of the predicted limit cycles can be obtained 
by substitution of the predicted limit cycle frequency into the 
"real part" given in equation (5-47) and solving the following 

Real Part (o)^ c ^ ' 

“ijjy 

Following this procedure the amplitude of the predicted limit 
cycles are shown in table 5-1. 


(5-45) 


(5-46) 


(5-47) 


(5-48) 
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Table 5-1. Limit Cycle Amplitude for Zero Telescope Look Angle 


LIMIT CYCLE 
FREQUENCY 
(Hz) 

RELATIVE RATE (!«„) 
LIMIT CYCLE AMPLITUD]! 
(rad/sec) 

RELATIVE ANGLE (i(Jy) 
LIMIT CYCLE AMPLITUDE 
(rad) 

5.206x10“^ 

5.393x10"^ 

1.649x10“^ 

1.419x10"^ 

1.048x10"^ 

1.176x10“^ 

0.136 

2.054x10"^ 

2.398x10“^ 

0.33 

4.809x10"^ 

2.317x10“^ 






One last piece of information to be determined is the stability 
of the limit cycle points that have been computed above. In order 
to ascertain limit cycle stability an analytical procedure outlined 
in a text titled ’’Multiple- Input Describing Functions and Nonlinear 
System Design” by A. Gelb and W. E. VanderVelde, is employed. This 
procedure requires that "jw” be substituted for ”S” in equation (5-42) 
and the resultant expression be written in terms of real and im- 
aginary parts as follows 

A = = 0 (5-49) 


In order for a stable limit cycle to exist the following condition 
must be met 






01 = 01 . 


2n 

doi 


L.S, 




L.S. 



(ll-W 



> 0 


(5-50) 


Performing the above outlined operations and computing the 
partial derivatives indicated in equation (5-50) gives 




=V r- 




[-4 


[9. 99371x10^0)^-9. 437xl0^wVl. 56xlO^(u^ 


.8795x10^0)^+1.4x10^0)^] 


] 




-6 (1. 5783xl0^°+9 . 99371xl0^N^^) 0)^+4 (2 . 2356x10^^+9 . 437xlO^N^^)o)^ 

(5-51) 


-2 (2. 87111x10^^+1. 56xl0^Nj^)o) 

^ *= -7(6.5039xl0®)o)^+5(7.283xl0^° +4.8795N^^y)a)^ 
-3 (3 . 57909x10^ Vl . 4xl0^Hj^y )U) ^2 . 6x10^ 
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Applying Uhe stability criterion of equation (5-50) shows that 

—3 

the limit cycles occurring at 5.206x10 and 0.136 Hz are unstable 

-2 

while those occurring at 1.419x10 and 0.33 Hz are stable* 

The actual limit cycles observed on the lOG time simulation for 
zero telescope look angle was at a frequency of 0*3229 llz and the 

amplitudes of the relative rate and angle were 6,375x10 ^ and 

3.142x10"*^ rad respectively* This Is in close agreement with the 
limit cycle frequency and amplitude predicted vici the analysis 
outlined above. However, a limit cycle at approximately 0*014 Hz 
that is predicted by the analy.sis was never observed on the actual 
time simulation. The reason for this is the finite bite size and 
sampling frequency that is employed in the simulation. The limit 
cycle amplitudes duo to the sampled data effects were approximately 
a factor 60 Icirgor than the level of limit cycle predicted at 0,014 
Hz by the anlaysis* Hence the system cannot sustain this low am- 
plitude limit cycle oscillation and goes into the limit cycle 
condition dictated by the sampled data effects. 

In addition, the analysis described above indicates that 
the limit cycle frequency is independent of both wire torque jump 
and vrire torque slope while ],imit cycle amplitude is directly 
proportional to wire torque jump. Both these effects have been 
verified on the lOG time simulation. 

5^1^3.2 Limit Cycles for 40^ Telescope Look Angle - Substituting 
the parameters listed in section 5.1.3 with = 40° and “ 

25 n-m/rad into equation (5-40) the following results. 

= 1.72907x10^5^3.96713x10^3-6.98984x10^“ 

5 

= 1. 30177x10^S^-C1.972x1oVN|^^)S-( 7. 0025x10^)- 

= 1.00613x10^5^-89.8955-187.5 

a, = 1.16036x10^*5^+1. 97 2x10^5+7 xloV 
4 h 

a” = 1. 03019 x 10^5 V(1.972xloVN^^y)S+(7. 0025 x10^)+ (5-52) 
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Eg = 3.4714x10^8^ 

= -2.91285x10^8^ 

a„ = 4.47752x10^8^-89.8958-187.5 

O 

= 2.878x10^8^+239.728+500 

Substituting equation (5-52) into equation (5-39) and evaluating yields 

A = 1. 91726xl0^^S®+3. 22055xl0^^S^+1.66498xl0^^s‘'+5. 05956x10^^8^+8. 46079x10^^8^ 
+8, 49024xl0^'^S^+2. 01736x10^ 27138x10^^8+1. 7x10^^+ — 

+N|^y [2 . 83152x10^^8^+1 . 4278x10^^^8^+3 . 0781x10^^8^+4 . 19x10^8^4 . 369xl0®8^ 
+1.1x10^8^+7.8x10^83 = 0 (5-53) 

Rearranging terms gives 

2 . 83152x10^^8^+1. 4278xlO^°sV3 . 0781xl0^°S^ 4.19x10^8^+4 . 369x10^8^+1.1x10^8^+7 . 8x10^8^ 

1. 91726x10^^8^+3.22055x10^^8^+1. 66498x10^^8^+5 . 05956x10^^8^+8, 46079x10^^8^ 

+8.49024x10^^8*^+2.01736x10^^8^1.27138x10^^8^1.7x10^^8+8.7x10^^ „ - ^ 

(5-54) 

Substituting s=jw into equation (5-54) results in 

+ 0 )^ [ (-2 . 83152xl0^°o)®+3 . 0781xlO^°to^-4 . 369x10^10^+7 . 8xl0^)+j (1. 4278x10^^(0^-4. 19xl0^t0^ 

+l.lxl0^to)] ^ 

(3 . 22055x10^^(0^-5 . 05956xl0^'*^(O^+8 . 49024x10^^(0^1 . 27138xl0"^ '^*(0^+8 . 7xl0^'*') 

+i (1. 91726x10^^(0^-1. 66498x10^ V+8. 46079x10^^(0— 2. 01736xlO^^(o\l. 7xl0^^o>) 
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(5-55) 


In order for equation (5-55) to be valid the imaginary part ot 
the loft hand of the equation iiuiat be equal to zero. This leads 
to the following constraint equation as a necessary condition for 
system limit cycles, 

w^^-79. 4585w^^+/*00. 298w^^-335. 99w*’+84 , 3163w^-7 .10739^^7 . 28899xlo"^w^ 
-6.79713x10“^ =0 (5-56) 


Solving equation (5-56) results in the following 
U) =* +.103175 
= +.032185 
= -,416312+j. 04158 
= +.416312+]. 04158 
= +2.0887 
= +.78875 
=+8.6092 


(5-57) 


In order to determine if each of the possible solutions shown 
in equation (5-51) are possible limit cycle points they each must be 
substituted into the real part of equation (5-55) and the result 
shown to be positive. The real part of equation (5-55) is given 
by equation (5-58) . 

Real Part = (-2. 83152x10^' V+3. 0781x10^ V‘-4. 369x10 V+7. 8x10^^ 

(3. 22055xl0^^u)®-5. 05956x10^^0)^^+8.49024x10^^0)^-1. 27138x10^^0)^+8. 7x10^^) 

10 5 9 3 7 12 9 14 7 

+(1.4278x10 0)^-4.19xl0 0) +1.1x10 0 )) (1.91726x10 0) -1.66498x10 0)^ 

+8 . 46079x10^^0)^-2 . 01736x10^^01^+1. 7xl0^^0)) 


' 2 . 

(3. 22055xlO^^O)®-5. 05956x10^^0)^+8. 49024x10^^*0)^-1. 27138x10^^0)^+8.7x10^^) 
+(1.91726xl0^'^o)^-1.66498xl0^^o)V8.46O79xl0^^o)^-2.G1736xl0^^o)^+1.7xl0^'^O))^ 


Substitution of the possible limit cycle frequencies shovm in 
equation (5-57) into equation (5-58) reveals that the solution for 
'V of 8,6092 rad/sec cannot be a limit cycle point. Following 
the same procedure as outlined in section 5. 1.3.1 the amplitude of 
the various possible limit cycles are given in table 5-2. 


The stability of the limit cycle points computed above is 
determined in the same manner as outlined in section 5. 1.3.1. The 
partial derivatives required to evaluate limit cycle stability 
are given in equation (5-59). 


2JL = -&[-2.83152xlO^°w®+3.0781xlO^*’a)^-4.369xloV+7.8xloV] 

t/lij 

y 



Vo 

■1) w 


■£)W 


[1 . 4278x10^ V-4 . 19xl0^w^+l. IxlO^w^ ] (5-59) 

’I' 

^y 

8 (3 . 22055x10^^+2 . 83152xlO^°N^y ) m^-6 (5 . 05956x10^ V 3 . 0781xl0^°N^y)a)^ 

+4 (8. 49024x10^^4 . 369xl0^K,^y )to^-2 (1. 27138x10^^7 . 8xl0^Nj^y)w 

9 (1. 91726xlO^^)to®-7 (1.66498x10^^+1. 4278xl0^°N^y)w^ 

I 

+5 (8.46079x10^^+ .4. 19x10 ^N^^y)ui^-3 (2. 01736xl0^Vl.lxl0^N^) w^+1. 7x10^^ 


Applying the limit cycle_stability criteria shows that the limit 
cycle points at 5.122x10 ^ and 0.126 Hz are unstable while the 
limit cycle points at 1.642x10”^ and 0.332 Hz are stable. 


The limit cycle observed on the lOG time simulation (see section 

5.2) was at 0.3226 Hz and the amplitudes of the relative rate (i|j ) and rela- 

_4 ~4 y 

tive angle (tJ^^) were 7.998x10 rad/sec and 3.946x10 rad respectively. 

The stable limit cycle at approximately 0.16 Hz predicted by the 
analysis was not observed on the time simulation for the reason 
outlined in section 5. 3.1.1. It can be seen from these results 
that the limit cycle frequencies arc hardly affected by telescope 
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Table 5-2, Limit Cycle Amplitude for 40“ Telescope Look Angle 


LIMIT CYCLE 
FREQUENCY 
(Hz) 

RELATIVE RATE (^y) 
LIMIT CYCLE AMPLITUDE 
(rad/sec) 

RELATIVE ANGLE (!()„) 
LIMIT CYCLE AMl’LITUDE 
(rad) 

5.122x10“^ 

5.901x10"^ 

1.833x10"^ 

-2 

“6 

-5 

1.642x10 

1.232x10 ° 

1.194x10 ^ 


-5 

-5 

0.126 

1.417x10 ^ 

1.797x10 


-4 


0.332 

6.693x10 ^ 

3.204x10 



look angle while the limit cycle amplitudes are slightly affected 
by the telescope look angle. In addition, if 1 arc-second total 
pointing stability is desired from the 106 and 20 percent of this 
budget is allocated for wire torques the jump in the wire torque 
characteristic should not exceed 0.05 n-m* For the small telescope 
described in section 4.1.4 the jump in the wire torque character- 
istic should not exceed 0*01 n-m if the same ground rules are 
followed. 


5.2 Nonlinear lOG Performance Characteristics - Figures 5r-3 

to figure 5-11 show 106 system time response for a 40 degree tele- 

scope look angle for a jump in wire torque characteristic = 1 n-m 

and a wire torque slope of “ 25 Figures 5-12 to 5-20 

Indicate lOG time response for a 40 degree telescope look angle with 
= 0.5 n-m and = 25 n-m/ rad. The legend used in these fig- 

ures is explained in table 4-1. Examination of these figures indicate 
that the telescope y axis limit cycle pointing error is directly pro- 
portional to the jump in the wire torque characteristic as implied 

in section 5.1* 3.2. The transients that are obsesrved at the start 

of the run are due to an initial wire torque of 17.45 n-m which is 

applied to the telescope when it is rotated to a forty degree look 
angle. Similarly it is seen that the y axis control torque limit 
cycles about this point since the control system must produce a 
torque equal and opposite to the average applied wire torque. In 
addition due to the geometry of the 106 and the fact that the tele- 
scope was only rotated about the y axis results in zero rotational 
motion of the telescope about the x and z axes, zero rotational 
motion of the pedestal about the x and z axes, zero translational 
motion of the pedestal CM along the y axis, and a control torque 
only being applied to the telescope y axis. 
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6* DERIVATION OF SUSPENSION CHARACTERISTICS FOR COINCIDENT 
GIMEAL SYSTEM CONFIGURATION 

The European Space Agency (ESA) in conjunction with Dornier 
Systems has recently modified the design of the lOG to the Coinci- 
dent Gimbal System configuration. The most significant change be- 
tween the lOG and the Coincident Gimbal is the configuration of 
the suspension system. The lOG had a square base pedestal with 
a four point suspension system consisting of a set of isolators 
at each corner of the pedestal* The pedestal for the Coincident 
Gimbal is circular with three sets of skewed isolators located 
on the periphery of the circle separated from each other by 120 
degrees. Each of the isolator sets consist of tv;o spring dampers 
skewed with respect to each other as shown in figure 6-1. In 
the sections that follow the coincident gimbal suspension linear 
and rotational stiffness and linear and rotational damping co- 
efficients are derived as a function of individual spring stiff- 
ness, damping, and soring slew angle. In the derivation it is 
assumed that all of the springs have equal stiffness and damp** 
ing characteristics along their respective axial directions. 

In addition, suspension system relationships inherent in 
the Coincident Gimbal design are derived and a comparison is made 
with the lOG four point suspension. 

6.1 Derivation of Linear and Rotational Stiffness and 
Damping Coefficients for the Coincident Gimbal Suspension System - 
The geometric configuration of the proposed Coincident Gimbal 
skewed spring suspension system is shown in figure 6'-l. 


Examination of figure 6-1 indicates that three spring damper assem- 
blies are positioned 120 cJcgrc=ies ap^irt around the periphery of a 
circle of radius R. Each assembly contains two springs that are 
canted at an angle a with respect to the plate (i.e*. Coincident Gimbal 
pedestal) to which they are attached. 


The force that is applied by the spring can be written 


as 


F . “ R[|r .+e[-lr * ' ] 

SJ * OJ ^ * OJ ’ 



where 


K - Axial spring stiffness (n/m) 
r . = Unstressed spring length 


(6-1) 


E = Spring elongation (i.e., translal Ion of pedestal C.M.) 
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Equation (6-1) can be written as 


F i == K{l(r .+e )^+(r ,+e +(r .+e )2]^^^_|;r2 +2 jl/2j 

sj oxj '• oyj ozj z'' oxj oy.j Oz.j ■’ 


[(r .+£; )^+(r ,+e )\(r ,+e )V^ 

OXJ x' ' oyj y' ozj z'' •' 


X {(r^ 4 +^„)i +(r •+£ )i +(r .-J-e )i 

OXJ x' X ' oyj y y ' ozj z' 


oyj y' y ozj 


(6-2) 


Equation ( 6 - 2 ) can be written as 


F . = K 
SJ 


1 -[r^ ,+r^ .+r^ .1^^^ 

OXJ oyj ozj^ 


I(r ,+e )^+(r ,+e )^+(r ,+t: )^]^^^ 

OXJ X oyj y^ • ozj z^ ^ 




Expanding equation ( 6 - 3 ) and eliminating second order terms 
(i.e., je|«l) results in 


F^.'«K{1- 7 -r ~7 r~\TT7:7H(r .+c )i +(r ,+c)i+(r ,+c )i } (6- 

sj 1 1+2 (r^^ ■ 1-^.^ . L^) ] 1/2 OXJ x^ x ^ oyj y" ozj y/ 


4\ 


^ 1^2 . 2 
r rrr .H’l 
OXJ oyj OZJ 


Making ui^o of the oxpnnslou of , vX/2 where x«l equation 

ti+x ) 


(6-4) can be written as 


F ,‘qc{r ,c 'h: ,c +r .e }{(r ,+G )1 +(r .+E )1 +(r ,+E )i } 

SJ ox.i X oyj y OZJ z ' oxj x x oyj y y ozj z z ( 6 - 5 ) 


2,2 ,2 

r .+r ,+r 

OXJ oyj OZJ 


F .= 


Eliminating second order terms from equation ( 6 - 5 ) gives 
K 


s j 2 .2 .2 

r ,+r .-IT 
OXJ oyj OZJ 


yi^ 

{(i' ■ ,r, +r .r .c +r .r .c )1 
OXJ >; OXJ oyj y oxj ozj z x 


+(r ,r .r +r ,e +r ,r .t, )1 
OXJ oyj X oyj y oyj ozj z 3 


+(r ,r ,c +r ,r ,e +r^ .e )i } 
OXJ OZJ X oyj OZJ y ozj z^ z 


( 6 — 6 ) 



For the 1,2 isolator pair sho\m in figure 6--1 the following applies 


r ,“~r sina 
osl o 

^oz2= 


The unstretched dimensions of Isolator pair 3,4 written in 
lOG pedestal coordinates can be obtained in the following manner 


^*0x3 


cos 120" 

sin 120" 

0 

^oy3 


-sin 120" 

cos 120" 

0 

_^oz3 


0 

0 

1 

^0x4 


" cos 120" 

+ain 120" 

0 

^oy4 


-sin 120" 

cos 120" 

0 

_^oz4 . 


0 

0 

1 


"oxl 


'oyl 


J 

- ^ozl J 


^^0x2 


^oy2 


_^oz2 _ 


Similarly the unstretclied dimensions of isolator pair 5,6 
in lOG pedestal coordinates can be written as 


^0x5 


“ cos 120" 

-sin 120" 

0 ' 


^0x2 

^oy5 


sin 120" 

cos 120" 

0 


’^oy2 

_^oz5 . 


0 

0 

1 


.’^oz2 _ 

^0x6 


' cos 120" 

-sin 120“ 

0 ~ 


’^oxl 



p 

0 





= 

sin 120 

cos 120 

0 


r_„. 


Substituting the appropriate dimensions for the unstretched 
spring length given for the various isolator pairs into equation 
(6-6) results in 

^s] (cos^a)c^i^+(sln^a)c^i^] 


( 6 - 10 ) 


F „ . = K{[-^os^ae + ]i +['^^^^os^ac + -^os^ac ]i ^'[2sia^ac } (6-11) 

s3,4 2 X 2 y X 2 x 2 y y a is 


F 


s5,6 X / y' X 

The total linear compliance force is given by 


^os^ac + |cos\ 

2 X 2 y' y 


= K{[-|cos^ar^~ -^^os2ac Jijl[ — ^^os'^ac^+ -^os^tL^Ji^ +[2sin^ae_^]i , } (6-12). . 
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F = F +F +F 
s si, 2 s3,4 x5,6 

F = K{[3cos^a]c i +[3cos^a]c i +[6sin^a]e i } 

s “■ ^xx yy •'zz 

Hence the linear spring constants using the same notation as defined 
in section 2*0 , 

4K = 4K = K = SKcos^a 
X y x,y 

4K = 6Ksin^a 
z 

Similarly, since the damping forces are applied in the same 
direction as the spring compliance force the fallowing applies 


(6-13) 


(6-14) 


4D == 4D ^ D = 3Dcos a 
X y K,y 

4D = 6Dsin^a 
z 


(6-15) 


where 

,D “ Damping coefficient along the axial spring direction ( — — — ) 

The torque applied by the suspension to the lOG pedestal due 
to pedestal rotation is given by 


T . - R, X F . 
S3 J S3 


(6-16) 


and 


e = 6 X R, 
3 


(6-17) 
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I 


Assuming has only an x and y component and substituting 

equation (6-17) into equation (6-6) and then substituting the results 
in equation (6-16) results in 


T . = - 


2 


X = -r — .R .0 -r R .0 +(r .r ,R .R ,-r .r .R .)0 ]i 

2 4^2 ,2 ozj yj X ozj xy yj y ^ oya ozj xj y:) oxj ozj yj z x 


r* 4 *t 4‘T* 

"cxj oyj ^ozj 


+[-r^ .R .R .0 +r^ .R^.6 +(r .r ,R .R ,-r ,r ,R^.)0 ]i 
OZJ X] yj X OZ3 xj y ' oxj ozj xj yj oyj ozj xj”^ z'^ y 

+[r ,r -R ,R ,-r ,r .R“.)0 +(r ,r ,R ,R ,-r ,r ,R^,)0 
oyj OZJ xj yj oxj ozj yj x '• oxj ozj xj yj oyj ozj xj'^ y 


+fr ,R ,-r ,R .) 0 ]i } 
oyj XJ oxj yj' z^ z 


(6-18) 


For the case under consideration 


R, ,, = -Ri 
1,2 y 

JiS* 1 

^3,4 = “ 

R. , = Ri + -I Rl 
5,6 2 X 2 y 


(6-19) 


Substituting equations (6—7), (6—10), (6-11), and (6-19) into equa- 
tion (6-18) gives 

^sl,2 “ 2R^K{(.,-in^a)0^i^+(cos^a)0^i^} C6-20) 

= r\{[ ^in^a0^+ ^sin^a0y]l^+ [| sin^a0y+ sin^a0^]iy+[2cos^a0^1i^} 


= R^K{[^in^a0^- sin^a0Ji^+ [| sin^a0„- -^In^a0^]i„+[2cos^a0,]i,} 




X y 


(6-22) 


Tlie total torque applied to the lOG pedestal is given by 

= T , .+T - ,+T . , 
s ,sl,2 s3,4 x5,6 


(6-23) 
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( 24 ) 






T = KR^{l3niiA]0^i +[3siii\4]0„l +l6cos^'a]e £ } 
s jv X y y £» fA 

Hence the rotationul stiffnens coefficients of the suspension 
arc given by 

k = k = 3iGl^sin^a 
X y 

k = 6KR^cos^a 

Z 

Similarly the rotational damping coefficients of the suspension are 
gi\ ?n by 

2 2 

d " d = 3 J:jR sin a 
X y 

2 2 

d = 6BR cos a 
z 


The parameters presently being proposed by Donier Systems for 
the suspension analyzed are listed below: 

K = 800 - 

III 


D = 


240 


n-m-sec 

rad 


a = 24.1 deg 
R - 0.25 m 

Substituting these parameters in equations 
(26) results in the following 

K = K = 2000 - 
X y tn 

K = 800 - 


z 


D = D 


in 


600 


n-ficc 


X y in 


( 14 ), 


(15), (25) and 


(25) 


(26) 
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D = 240 
z m 


k « k « 25 ^ 
X y rad 

k = 250 
2 rad 


d = d =7.5 
X y 


n-m-sec 


rad 


j _ - 7 C n-m-sec 

d - 75 1 — 

z rad 


6.2 Inherent Relationships Between Suspension Parameters for 
the Coincident Gimbal System 


Examination of equations (6—14), (6—15), (6—20) and (6—26) indicate that 
there are only four independent parameters that can be used to set all 
of the suspension characteristics. Since there are twelve parameters 
(i.e., linear stiffness and damping along three axes, rotational stiff- 
ness and damping about three axes) to be set for the suspension it is 
clear that these cannot be set independently with only four parameters. 

In addition, the linear stiffness along the three axes cannot be set 
independently since the stiffness along the x and y axes must always 
be equal with the stiffness along z set to some desired value by ad- 
justment of the angle “a”. The linear damping along the x, y, and z 
axes can only be controlled by varying the damping coefficient 
Examination of equation (6-15) indicates that the damping coefficient 
along the x and y axes are always equal but can be set to a desired 
level. However, the damping coefficient along the z axis is now de- 
termined and cannot be controlled. 


Assuming that the damping coefficient ”D" is set so that a 
particular damping ratio "c" is obtained about the x and y axes, 
this can be written as 


"x,y 


1 

2 




K 


m 


(6-27) 


The damping ratio along the z axis is given by 




z 


1 

2 


D 

z 



(6-28) 




> 



However 

D = 2tan^aD (6-29) 

K = 2fan^aK 
z x,y 

Substituting equation (6-29) into equation (6-28) results in 
C = a/ 2 tanaCo (6-30) 

Examination of equation (6-30) Indicates that the linear damp- 
ing ratio along the z-axis will equal that of t!ie x and y axes 
for a = 35.26 deg. For va]ucs of ”x" below tliis value the damping 
along the z-axis decreases \i?hilc it increases for cx above 35.26 deg. 

If the a angle is set at 15 deg. or bel iw, appreciable reduction 
in linear damping along the z-axis results which can result in 
deteriorated system response characteristics. 

The rotational stiffness of the mount about the x and y axes 
can be set by varying the distance "R” as indicated in equation 
(6-25). However, the rotational stiffness about the "z" axis is 
determined and cannot be controlled. In addition, the rotational 
damping coefficients about the x, y, and z axes are now determined 
and cannot be separately controlled. 


It will be shown (section 7) that the critical suspension parameters 
that govern lOG pointing performance are the rotational stiffness about the 
X and y axes and the translational stiffness along the z axis. Hence, 
one might conclude that increasing the linear stiffness along the x and 
y axes would result in a redact Lon of the radius of the base of the 
Coincident Cimbal pedestal since it seems reasonable that tlie J.ongi— 
tudinal stiffness K of the isol/itors would increase. However, if equa- 
tions (6—14) are divided by each other the foilowing results 


tana = 


K 


'2K 


x,y 


1/2 


(6-31) 


Equation (6-31) implies the geometrical relationship shown in figure 6-2, 


The value of K from equation (6-14) and the geometrical relationships 
shown in figure 6-2 can be written as 

K 


K = 


K 

x.y ^ 
2 

3cos a 


JhX. 


2K 

. 

+21C ' 

z x,y 


2K +K 

x.y z 


(6-32) 




6-8 


Substituting equation (6-32) and the geometrical relationships into 
equation (6-25) yields: 


k=k=k =3(- 
X y x,y 


2K +K 

x,y a 


K 


31 - 


2K +K 
x,y z 


2 2 
]R^ = (— )R 


(6-33) 


therefore 


R 


2k 

'■ K ^ 


(6-34) 


Equation (6-34) indicates that the radius of the base of the 
Coincident Gimbal ped€'Stal is completely determined by the rotational 
stiffness about the x and y axes and the linear stiffness along the 
z axis. Therefore, varying the linear stiffness along the x and y 
axes will not affect the radius of the pedestal base and cannot be 
used to affect a size reduction. 


6.3 Comparison of Coincident Gimbal Suspension System with 
Square Based Suspension System - A comparison can be made as to the 
required area of the pedestal base for the Coincident Gimbal Suspension 
System and the corresponding square base suspension. Assuming that 
tiie rotational stiffness about the x and y axes and the translational 
stiffness along the z axis are equivalent for both systems, the side 
of the square base suspension system can be written as 


2a 


120 


where 


k 

2(^) 


1/2 


(6-35) 


2 ai 2 o ~ length of one side of the square base 


The area required by the square base suspension is then given by 


4k 


S.B. 




K 


(6-36) 


The area of the corresponding Coincident Gimbal pedestal base 
using equation (6-34) is given by 


2k 


A. . = Tr(-r^) 


C.G. 


or the ratio of the area of the Coincident Gimbal Base to the square 
base is 


6-9 


~ - 1.571 (6-37) 

Hence, the area of the Coincident Gimbal pedestal base is almost 60 
percent greater than the required area for the square base suspension 
system. 


A. 


C.G. 


S.B. 







7. COINCIDENT GIMBAL SYSTEM PERFOEMANCE 

In the sections that follow, the Coincident Glmbal system per- 
formance is described for the nominal suspension parameters speci- 
fied by Dornier Systems. Both stellar pointing and earth tracking 
performance are discussed for the given suspension system. Based 
on the results obtained, the suspension parameters are modified in 
order to improve the performance of the Coincident Gimbal System 
in both stellar pointing and earth tracking. Representative 
system performance characteristics are shown for the modified 
suspension system and a recommendation is made as to the suspension 
parameters that should be implemented for the Coincident Gimbal 
System. 


7.1 Coincident Gimbal Pointing Performance - In order to ob 
tain the Coincident Gimbal geometrical system configuration the 
following parameters must be modified from the values shorn in 
table 2-1 as outlined below: 

= 138.45 Kg 

■^2x "" 12.695 Kg-m^ 

= 16.313 Kg-m^ 

J 22 = 5.417 Kg-m^ 

Rg 2 o = .1017 i^ - .0029 iy - .5251 i^ m 

= -4.822 i + 0.0029 i - 1.429 i m 
120 X y . z 

r, = 0.1017 i - 0.0029 i + .275 i m 
1 X y z 

System time responses were taken for the suspension para- 
meters furnished by Domier which entail changing the baseline 
parameters in the following manner. 


4D = 

4D 

= 600 

X 

y 

m 

4D = 

240 

n-sec 

z 


m 

11 

4K 

y 

= 2000 n/m 

li 

800 

n/m 
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d 




d 


k 


X 


- d =7*5 n--m-"Sec 

y 

"75 ri"m*-sec 



25 


n-m 

sec 


k 

z 


^ 250 


a-m 

rad 


For a telescope look angle of sixty degrees and a crew motion 
disturbance applied along the z axis the peak pointing error in- 
curred about the telescope y axis was 4#34x]0“6 rad (.895 arc- 
second) . This value of pointing error from one source is too 
large if an arc-second of pointing stability is to be realized 
when considering all of the other contributing system error 
sources (e,g., sensor and actuator noise, gimbal nonlinearities, 
sampling and quantization, etc*). Hence, it is clear that the 
y axis pointing error must be reduced. Since it is not advisable 
to increase the pointing control loop bandwidth beyond 2 Hz from 
system noise and structural flexibility considerations the only 
way to reduce this error is to reduce the suspension linear stiff- 
ness in the z direction to at least the baseline set of parameters 
which would yield similar pointing performance as that obtained 
for the lOG as outlined in section 4. 

Figures 7-1 through 7-21 show system pointing performance 
for a telescope look angle of zero degrees with the crew motion 
disturbance shown in figure 4-1 applied simultaneously along the 
X, y and z axes of the shuttle orbiter for the Coincident Gimbal 
System parameters specified by Hornier Systems. Figures 7-22 
through 7-42 and figures 7-43 through 7-63 show Coincident Gimbal 
pointing performance for telescope look angles of 40 and 60 degrees 
respectively with all other parameters being the same as above. 

The legend used in these figures are explained in table 4-1. 

7.2 Coin::ident Gimbal Earth Tracking Performance - Using 
the earth tracking command profiles described in section 4*2, the 
perfomance of the Coincident Gimbal System as defined by Hornier 
Systems was determined* The result of this determination are shown 
in figures 7-64 through 7-84. A definition of the nomenclature 
used in these plots are given in table 4-1. Examination of figure 
7-66 indicates that the peak tracking error incurred was 0.262x10"^ 
rad (.540 arc-second) while an examination of figures 7-72, 7-77 
and 7-79 indicate that the pedestal y rotation was 0.311 rad 
(17.8 deg) with the x and y translation of the pedestal C.M* equal 
to 0.164 m (6.46 in) and .034 m (1*34 in) respectively. These 
values correspond to an isolator elongation of 4.169x10“*^ m (1.641 in) 
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which is too large and could not be accommodated within the frame- 
work of a reasonable design. This situation becomes worse as 
altitude of the orbit is reduced from 438 KM for which the above 
values apply. Hence it is clear that if the IPS is to have a reason- 
able earth tracking capability the suspension parameters proposed 
by Hornier should be modified by stiffening the suspension rota- 
tionally and translationally as much as possible while still meet- 
ing stellar pointing requirements. This approach will be fully 
discussed in section 7.4. 

7.3 Coincident Gimbal Raster Scanning Performance - Examina- 
tion of experiments proposed to be integrated with the IPS indicate 
that many have raster scan requirements. In order to examine^ in 
a preliminary fashion, the raster scan capability of the Coinci- 
dent Gimbal System the raster scan performance of a particular 
experiment was determined. The experiment chosen for this 
evaluation was the X-UV (soft x-ray) Spectroheliograph. A schem- 
atic diagram of this telescope with its associated mass character- 
istics is shown in figure 7-85. 

The raster scan field that was considered was 90 x 90 arc- 
second with a scan time per line of 1.12 arc-second which yields 
a scan rate of 3.897x10“^ rad/sec (80.36 arc-second/ second) . The 
desired accuracy was that the tracking error should be within 
+1.7 arc-second for ninety percent of the time it takes to scan 
a line (i.e., 1,12 second). The technique used to accomplish the 
raster scan profile was to use simultaneous rate and position 
commands to the Coincident Gimbal pointing control system. The 
rate and position command profiles employed are depicted in fig- 
ure 7—86. 

Two computer runs were made for the X-UV Spectroheliograph 
using the above profiles. The first was to raster scan the tele- 
scope about the y rxis only. The second was to raster scan the 
telescope about an axis in the x,y plane making inclined 45 de- 
grees with respect to the x and y coordinate axes. For both 
cases the telescope was positioned so as to allow a symmetrical 
scan about zero look angle, and the pointing control loop band- 
width was set at 2 Hz. 

The results of these studies indicate that the system did 
not meet the accuracy specification of 1.7 arc-second for ninety 
percent of the line scan time for both cases Investigated. The 
tracking errors remained above 1.7 arc-second for approximately 
0.82 second of the 1.12 second that it takes to scan one line 
with a peak tracking error of approximately 9 arc-seconds. 

However, the peak control torque reqviired was 2.89 n-m which 
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is well within the torque capability of the Coincident Gimbal 
torque motor capability. These results indicate that raster 
scanning to the accuracy described above cannot be met by the 
"simple” scan profiles employed. It is anticipated, however, 
that if a torqu'^ command used in conjunction with consistent rate 
and position command profiles were employed, system raster scan- 
ning performance would be greatly improved and the accuracies 
desired would probably be achieved* It is recommended that future 
effort should be expended to implement this raster scan command 
technique and determine the accuracies that can be achieved. 

7.4 Pointing Performance with Modified Suspension Character- 
istics “ It has been shown in sections 7.1 and 7.2 that "large" 
pointing errors (i.e., 0.895 arc-second) result in the presence 
of crew motion disturbances while excessive pedestal rotations and 
translations with accompanying "large" (i.e., 1.7 in) isolator 
elongations occur during earth tracking translations for the iso- 
lator parameters defined by Dornler Systems, Hence, it is evident 
that the suspension parameters should be modified in order to give 
the IPS a reasonable slewing capability while still meeting stellar 
pointing performance requirements. It has been demonstrated that 
increasing the rotational stiffness of the suspension system about 
the X and y axes reduces the pedestal rotations and translation 
from those incurred using the nominal suspension parameters for 
the lOG system. Hence, it would similarly be desirable to increase 
the rotational stiffness of the coincident gimbal suspension over 
those specified by Dornier xijhich are approximately the same as 
those of nominal lOG in order Co minimize pedesral rotations and 
translations incurred during slewing. However, che maximum radius 
Chat the coincident gimbal base can have from mounting considera- 
tions is 0.9 meters. Since from stellar pointing considerations 
it would be advisable to reduce the translational stiffness along 
the z axis to 500 n/m in order to reduce the pointing error in- 
curred due to crex^ motion disturbances, the maximum rotational 
stiffness that could be realized about the x and y axes is 202.5 
n-m/rad (see section 6,2)- 

The other suspension parameter that can be varied is the 
translational stiffness along the x and y axes keeping the trans- 
lational stiffness along the z axis equal to 500 n/m and the rota- 
tional stiffness about the x and y axes equal to 202.5 n^rm/rad. 

The effect of this parameterization on peak pointing error for a 
telescope look angle of 60 degrees is outlined in figures 7-87 and 
7-88. It should be noted that in the parameterization of the linear 
suspension stiffness along the x and y axes the linear suspension 
damping ratio was maintained at 0,1 along these axes when consider- 
ing all of the mass above the suspension point (i.e., m^i-m^). 
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Examination of these figures indicates that the y axis pointing 
error decreases from 0.866 arc-second to 0.555 arc-second when 
the linear suspension stiffness is increased from 500 n/m to 
10,000 n/m along the x and y axes for a crew motion disturbance 
applied along the z axis. Examination of figure 7-25 indicates 
that the peak pointing error incurred about the x axis is inde- 
pendent of the linear stiffness along the x and y axes further 
demonstrating that the suspension system does not act primarily 
as a classical isolator for these errors as described in section 
4.1. It should also be noted that most of the pointing improve- 
ment is realized when the suspension linear stiffness along the 
X and y axes is Increased to between 2000 and 4000 n/m. Hence, 
in order to minimize the changes to the suspension parameters 
specified by Domier while still realizing most of the advantages 
that can be achieved within the constraints of the coincident 
glmbal system the recommended suspension parameters are listed 
in table 7-1. 

The suspension parameter values given in table 7-1 result in 
the following Coincident Glmbal suspension system isolator design 
characteristics . 

K = 750 n/m 

D = 178.2 n-sec/m 

a = 19.47 deg 

R = 0.9 m 

System time responses for the recommended suspension para- 
meters with crew motion disturbances applied along the x, y and 
z shuttle orbiter axes simultaneously are given in figures 7-89 
through 7-109, 7-110 through 7-130, and 7-131 through 7-151 for 
telescope look angles of zero, 40, and 60 degrees respectively. • 

The nomenclature used in these plots is explained in table 4-1. 
Examination of the plots for the angular rotation and the trans- 
lation of the CM of the pedestal Indicate that they do not seem to 
be returning to zero. This is in fact the case since the crew motion 
disturbance, which is mmomentum conservative, will still cause the 
shuttle to rotate through some small angle from its initial start- 
ing point. This will then cause the pedestal to rotate through 
the same angle once the transients die down. Since the transla- 
tion hhe pedestal C.M. is measured with respect to Rj^20 

which is an inertial vector it will also approach a nonzero value 
once transients die down due to shuttle orbiter rotation. 
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7.5 Coincident Gimbal Earth Tracking Perforroance with Re- 
commended Suspension Parameters Using the earth tracking com- 
mand profiles described in section 4.2 the performance of Coinci- 
dent Giiubal system with the recommended suspension parameters 
was determined. These results are shown in figures 7-152 through 
7-172. A definition of the nomenclature used in these figures 
is given in table 4-1. Examination of figures 7-90 indicates that 
the peak tracking error incurred was 0. 352xl0“^rad (.726 arc-second) 
which is withint the desired +1 arc-second tracking accuracy but is 
approximately 34 percent worse than that realized for the suspen- 
sion parameter specified by Dornier. However, the peak pedestal 
rotation was 0.0273 rad (1.564 deg) and the peak pedestal C.M. 
translations were 15-6 mm (.61 in) and 2.42 mm (.095 in) along 
the X and z axes respectively. These values are more than an order 
of magnitude smaller than those incurred for the suspension para- 
meters defined by Dornier. These values will result in isolator 
elongations of approximately 8.5 ram (.0335 in) which can easily 
be accommodated by the Coincident Gimbal suspension system design* 
Hence, it has been demonstrated that the recommended suspension 
parameters not only meet the stellar pointing stability require- 
ments of + 1 arc-second in the presence of crew motion disturb- 
ances, but also greatly reduces the pedestal rotations and trans- 
lations and hence isolator elongations incurred during earth 
tracking making these values consistent with a reasonable sus- 
pension design. 

Additionally, a comparison between system response for the 
recommended Dornier suspension parameters (figures 7-65 through 
7-85) and that obtained for the recommended suspension parameters 
indicate that the system responses are considerably different in 
functional form. The difference in response is due to the stiffen- 
ing of the suspension system for the recommended parameters which 
results in the more oscillatory response characteristics for this 
system. However, the increased oscillatory nature of the response 
has no detrimental effect on overall system performance. 
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Figure 7-1. Telescope x-Axls PolnClng Error 
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Figure 7-5. Telescope y-Axls Rate Error 
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Figure 7-6. Telescope z-Axis Rate Error 
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Figure 7-fl. Pedestal y-Axia Rotation 
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Figure 7-9. Pedestal z-Axls Rotation 
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Figure 7-11. Pedescal y-Axla Rate 
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Figure 7-12. Pedestal z-Axis Rate 
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Figure 7-13. Pedestal x-Axls CM Translation 
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Figure 7-17. 
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Figure 7-23* Telescope y-Axls Pointing Error 
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Figure 7-25. Telescope x-Axls Rate Error 
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TELESCOPE LOOK ANCLE - 40* 

CREW MOTION Dl. ;URBANCE ALONG THREE AXES 


TIME (SEC) 

Figure 7-30. Pedestal z-Axis Rotation 
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Figure 7-31. Pedestal x-Axls Rate 
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Figure 7-35. Pedestal y-Axls CM Translation 
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Figure 7-37. Pedestal x-Axia CM Velocity 
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Figure 7-38. Pedeatel y-Axia CM Velocity 
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Figure 7-41. y-Axia Concroi Torque 









THTYST <.RAD)*10-^ THTXST <.RAD)*1 

-0.02 0.00 0.02 0.04 0.06 -0.20 -0.10 0.00 0.10 i 



TELESCOPE LOOK ANCLE - 60* 

CREW NOTION DISTURBANCES ALONG THREE AXES 


Figure 7-A3. Telescope x-Axis Pointing Error 
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Figure 7-44. Telescope y-Axls Pointing Error 
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Figure 7-46. Telescope x-Axie Rate Error 
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Figure 7-67. Telescope y-Axls Race Error 


TELESCOPE r.OUK ANCLE - 60 * 

CREW MOTION DISTURBANCES ALONG TMRE;E AXES 
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Figure 7-56. Pedestal y-Axls CM Translation 
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Figure 7-59. Pedestal y-Axla CM Velocity 
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Figure 7-60. Pedestal *-Axls CM Velocity 
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Figure 7-61. x-Axls Control Torque 
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Figure 7-62. y-Axia Control Torque 
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Figure 7**68. Telescope y-Axls Treclclng Rate Error 
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Figure 7-79. Pedestal x-Axls CM Velocity 
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Figure 7-80. Pedestal y-Axie CM Velocity 
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Figure 7-81. Pedeatal r-Axie CM Velocity 
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Figure 7-83. y-Axia Control Torque 















Figure 7-85. X-UV (Soft X-Ray) Spectroheliograph 
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Figure 7-88. Telescope Peak Pointing Error vs Linear Suspension Stiffness 


r %t!MA ir)M«| tw ril twHI liowt fen dii l ImM 


THTYST CRADJ^IJSr-s THTXST < RAD >^10-3 







RAD/SEC )3Sl0-5 





OMSZST < RAD/SEC) a-: 1J0T-13 OfJGVST CRAB/ SEC 






u«iK MMiu: - ci" 

cm'wni'i'inn iiiimuiiiANCi-; Ai.nut: au. iittiiu- asiijj 

itKCfii’U'U Nni.ii conurini’Mi* cimkai. rAuAiti mis 



I Cxrm 'ci.ui 'ip.a-f 

TiMf 

FAgutrc: 7-^5* Pudciital x-Axls KiilnLjuii 


FVLIN 


TKI,KS(‘«n'K unciK ANCl.F - (l“ 

CHEW Hm'inti niSTnitHAKt:E Ai.tmt; all thiii-.e axes 

RKCUHMLNUHIJ COtHElDENT iNMIlAL .StiSl'I-JNS LnH r'AHAM^n’E!lS 




OMZGABX <;RAD)5K10-2 THETA2Z <;RAD)5Klj 






c.,nrr 

7“^*^ 


nixu tin anj;u- ^ n” 

(ittr.w H": tns’ii’i!HAr;i;i: auinc ai.u iiim-K axi’.s 

CMiK'ijni-ur cihhak I'AHAhi - j 



!:I,KI llt.tUi 'irMUf \a,W Vi.KI 

TH'H. f XliC ) iniit f 1 -I 

y-Ai(iH Hid 0 


TELKSCOI'K !,«OK AN(ILK “ 0* 

CJtEW MOTKW DJS'ruiIjlANCR AlAlNC ALL THHl-E AXKS 
REtlOMHKN|)i:|) COJtICIliUNT CIHKAL SUill’i-NliUlK FAKAMFTKKS 




EPS12Y < METERS >5)^lJEr-e EPS12X < METERS 






0.24 







RUN 213 












ORIGUij'^L' PAGE IS 
OE POOE QUALOT 


7-65 











RUN PM 


Figure 7-1]^. 


TIME CSEO 

Telescope y-Axis Rate Error 


TKLKSCDEE UIOK AN(JEK « AO® 

CHKW MOTltlN mSTUKRANCK AUlNi; TIIKKE AXES 
TilvCOUMKNDKlJ SUSPENSION PAliAMETKRS 








RUN 21*1 


Figuta 


TIME (SEC) 

^-118. Pedestal z-Aatls Rocatlon 


'iKiKScni^t: inni: Amiii-; “ 

CUKW MIITItm DISttiKllAHCfi AU^Hl! TllUKK AXKS 
RKCOMMKtmKM SUill'KNKKm fAKAHUTKIlS 








0P1EGA2Z <.RftD)?HlJ2r-3 0MEGA2Y CRAD^s^l^-^ 












a0“3 





















Figure 7"“136 


TIME <SEC) 

Telescope y-Axls Rate Error 


RUN 215 


TELESCOPE LOOK ANCLE - 60* 

CREW KCITION DISTURiiANCE ALONG THKEE AXES 
RECOMMtliUEn SUSPENSION PARAHICTEltS 



8*js0T iGf.m la.m i6,w 

TIME (SEC) RUN 215 

Telescope z-Axis Rate Error 



7-79 




Tint ':se:.c> 

Pedestal x-Axis Kotatlon 


RUN 51*5 


TKLliSCOPU I.tlOK ANCI.K - fiO" 

CRi:W HDTKIM Ji I STUUllAIiCK AUIJICJ TMRKK AXKf: 
RECUMKI-IMUK]) SllSPi;«SIU(f 1 'ARAMfi'J‘KHS 







0C1EGA2X ';RAB)3KlJBf-2 THETA22 '^RAD)ssi0“3 












page is 



EPS12Y (METERS EPS12X ( METERS 0-2 










EPS12ZDT ( METERS/SEC EPS12YDT ( METERS/SEC )5*fl0~3 






THTYC < NEUTON-flETERS ) 

,60 -0.<i0 -e.a? je .00 0.2£r 0.60 


6.00 8.00' 10.00 

TIME <SEC > 

Figure 7-149. x-Axle Control Torque 


RUN sn 






THTZC (HEUTON“NETERS)*I0-8 





YST 








m-ctmmiMDiii) susi»khri<in pajiamktkus 
KA irni vom TUACKimi h'wm r» 3 B km oiibit 

UBimi APPHOKUIATUJN TO KAHTIJ TRACKING PROFILES 








(VV-'A. 






11j*qh 6Qr.f;a^ 75. ej W.Cd^r 
TIME <SEC) 

Figure 7-157 •; telescope z-Axls Tracking Kate Error 


Wo.OUf 


K?Q.J 3 tf 


RUN 


T^,j00t m "etr.jaiar ys^so* % 

time: <se:c) 

Figure 7 - 156 . Telescope y-Aitls Tracking Rata Error 


sff.m 


10 ^ 5 , lakT 


W.tiEr 


RUN dFA 


RKCDMNl'tniKn SUSPKNSION l‘AllAHh:TKHS 

KAItTl! I^IHHT TJlACKlfJO nUlM JH KM OltiUT 

liSIWG AIM'ROXIHATJDN TO EARTH TRACK INC PUOFl LKS 


* M 


*1 «l w ir.lf/M VM'jV W 4* 


THETA2V XRAB THETA2X CRAP >5^10-^ 






RECOMMENDED SUSPENSION PARAMETERS 




3Qr.B0T 6S.m 75*00' 

TIME (SEC) 

Figare 7-159* Pedestal y-Axle Rotatlen 


90*00 105.00 120.00 

RUM 221 


7-91 











0MEGA2H (.RAD >5^10-5 0MEGA2Y <.RAD)s«10*i 


RECOHHE«DED SUSPEKSIOH PAUAJ1KTER5 

EARTH POINT TRACKING FROM A3fl KH ORUIT 

USING APPROXIHATIOM TO EARTH TRACKING PROFILES 


.m 45. £H W.m 75.J03* 

TIME CSEC) 

Figure 7-162. Pedeacal y-Axis Rote 


105*BQT isor.iaar 
RUN 2El 


HECOHMENJmi) RlJSPENSIClN PARAHKTERS 

EARTH POINT TUACKINC; FROM A38 KM OHBIT 

USING APPUOXlHATroN TO EARTH TRACKING PROFILES 




^5.J3J? W,m 71 

5.J33* 

9i3‘.00' 

iss.tjg 

iE0‘.iaar 

TIME CSEC) 



RUN 

221 

Figure 7-163. Pedpstol r-Axla 

Rote 









EPSISXDT < METERS/SEC 0-2 EPS12Z C METERS 

-S'. 16 -0.B3 S.B0' S.S3 0.16 0.E4 0.3E ^,16 IE -0.08 -0.04 S.SST 



7-95 




EPSISZDT (METERS/ SEC >H-clJBr-2 EPSIEYDT C METERS/SEC)5 k10-5 







THTYC < NEWTON-METERS) „ THTXC t NEWTON-METERS )M5lJ0r-2 




THTZC < NEWTON-METERS 

-e.asr -0,10 0.00 0.10 0.20 .1 


RECCMMiiNDRD SUSPENSION PARAMETERS 

EARTN POINT TRACKING FROM 43S KH ORUIT 

USING AI'PROXIKATION TO EARTH TRACKING PROFILES 


OJ 


i. 




B.m 


\5,ja0f 'w.w 




• J 3 [,m 75. te 

tine; (se:c) 

Figure 7-172 4 z-Aicis Control Torque 






RUN 231 


I- 

t 


r 

k 


:i=: 





8. ANALYSIS OF EARTH TRACKING PROFILES 


In this section the profiles that are required to perform 
earth point tracking from earth orbit will be derived. In this 
derivation a circular orbit and spherical earth are assumed and 
the effects of earth rotation is neglected. The maximum accelera- 
tion points for the earth tracking profile will be determined as 
a function of the various orbital parameters and an approximate 
formula will be derived which gives peak acceleration as a func- 
tion of orbital altitude. In addition, the peak rate of change 
of acceleration will be determined as a function of orbital 
parameters from which an approximate relation will be derived 
which gives peak rate of change of acceleration as a function 
of orbital altitude. These approximate relationships will be 
used to explain the behavior of the Coincident Gimbal System 
earth tracking performance characteristics, given in section 9.0, 

8.1 Derivation of Earth Tracking Profiles - Figure 8-1 is 
a schematic diagram of a vehicle in a circular earth orbit. 

The coordinate frame shown in figure 8-1 is inertial with 
the y axis along the orbit normal and the x and z axes in the 
orbital plane. Since the effects of earth rotation are being 
neglected in what follows, it can be assumed that the inertial 
coordinate frame is attached to the earth with its origin located 
at the earth’s center. Using the vector quantities shown in 
figure 8-1 a uiiit vector along the line of sight from the vehicle 
to some arbitrary point on earth can be written as 


A . 7 _ 

lf?l ~ 


( 2 - 1 ) 


where 

i g = Unit vector along the line of sight from the vehicle to an 
* * arbitrary point on the earth surface, 

- Vector from center of earth to point on earth's surface 

R - Vector from center of earth to vehicle, 
o 

Assuming a constant orbital rate about the positive y axis 
can be expressed as 


y ^ components of expressed geocentric 

coordinates 


R “ Radius of the earth 
e 

h = Orbital altitude 

0 ) = Orbital Rate 

o 

Expressing IL as 




and substituting equations (8-2) and (8-4) into (8-1) gives 


" (Pg ^ h]co5 QJ ?x ^ Vl ^ (^e S^nOjl^ 


Letting 

(8-b) 
fa-7) 

The inertial rata of change of the unit vector along the 
line of sight can be written as 


A = Xl " (/?c + M cos 6 

B “ Zj -f (Rc-fb) SfnB 
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■ji 


[ ( AB - A b) B Yt* a] - Yt f A A + B a) + [ Xf 8 ^ (A8 - A 8) a] 


(A*+Y’ + e®) 


V2 


( 8 - 8 ) 


The inertial rate of change of the unit vector along the 
line of sight from the vehicle to a point on the earth can also 
be written as 


where 


di 
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lx 


^ "1 







tOy 
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(8-9) 


to = An inertial rate written in geocentric coordinates. 
Expanding equation (8-9) gives 
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^'iSx 


- CO 


H iiSa) ty 


+ f ^iSy - A)y i^sx ) (8-/o) 


Comparing terms in equation (8-10) with those in equation C8-8) 
yields the following for the inertial rates required for earth 
point tracking written in geocentric coordinates 


to, = 


YB 


A'^4- Y^^+ B = 


{8-jO 


tOy 


(AB-AB) 
A' ■* y/ + 


( 8 -/ 2 ) 
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-Y,A 

A"' * Y,.^ + 8^ 


(8-/3) 




and 


A = 

^ fl) N f B 

(8-f4) 

8 = 

05^ 

(9-/S-) 


It should be noted that if the body rates required to perform 
earth tracking would be desired for an arbitrary telescope orienta- 
tion the rates given in equations (8-11) through (8-13) would have 
to be transferred to telescope coordinates and the transformation 
would need to be continually updated via a strapdown equation 
formulation. 

Assuming the point on earth to be tracked is contained in 
the orbital plane (which would yield maximum tracking rates and 
accelerations) and assuming that the earth point is along the 
X axis of the geocentric coordinate frame shown in figure 8-1 
gives the following relationships 

(8-U,) 

= (8-/7) 

Substituting equations (8-6), (8-7), (8-14), (8-15), (8-16) 
and (8-17) into equations (8—11) through (8—13) yields 


w, -- ~ ^ 


1 8 - /8) 




fffie - ffp f Rv * />) TrvS B} 




( 8 - 18 ) 


The line of sight angle about the y axis can be obtained by 
integration of equation (8—19) or by derivation with the aid of 
figure 8-2. 

Referring to figure 8-2, the following relationship applies 
between the variables 


/?.4h ^ 

Cin[l 8 ')-(i 80 - 


8-4 


sir, (iBO-ip'j 


(8-2^') 


I ) 


I ^ 

t i 

i : 


or 


Re + h 

Sin fp 


Sfn 


(p-e) 


(e-2i) 


where 

(J) = Line of sight angle and is equal to the direct integral 
of 0) 

y 

Solving equation (8—21) for ^ yields 


UfieU 




l^)cos 0 ^ /?e J 


(8-22) 


Equations (8-19) and (8-22) are the functional form of the 
angular rate and position command profiles (i*a*j and 

that are employed to command the telescope attitude control system 
to accomplish earth tracking. These profiles must be modified by 
the appropriate constants to reflect the desired telescope slew 
axis (i«e *5 they get multiplied by the components of a unit vector 
along the desired slew axis) which enables telescope slewing about 
any axis or alternatively allows the telescope to be in any orienta- 
tion mth respect to the orbit ihormal. 

The angular acceleration required to earth point track is 
obtained by differentiating equation (8-19) with the folloi-Ting 
result 


f ■ 




P„h(Rci-b)(2R^-i‘}^) sin 9 


-I 


fi) “ 2 ^?e coi 0 J 




( 8 - 2 ?) 


To obtain the points at which peak acceleration occurs, 
equation (8-23) is differentiated with respect to time and the 
results set equal to zero which gives the following 


Wy = -iOoPeh(^«4M(2/?e+>.)| J 


(8-24) 


8-5 




Setting equation (8~24) equal to zero results in the following 

2Re(fie+OcOS*e + [R« + O?«+h)*]cos0 (8-2S-) 

Solving equation (8-25) for cos0 gives the following consistent 
solution for the peak acceleration points. 


-[/?e + +32/?e 

COS0S 

h) 


( S - 2 6) 


Notice that two solutions exist for 0, one positive, the 
other negative which will yield two equal and opposite peak ac- 
celeration points equally positioned about the point where the 
vehicle is directly overhead (i.e., 0-0). 

The equations given above for the tracking command angular 

position, rates, accelerations and rate of change of acceleration 

involve the orbital rate w which in itself is a function of 

o 

orbital altitude. The relationship between (fl and orbital alti- 

o 

tude can be obtained in the following manner. 

For an orbital vehicle in a circular orbit assuming a spheri- 
cal earth, the following relationship applies 




(8-27) 


where 


G “ Gravitational Constant 

m = Mass of the earth 
e 

m = Mass of the vehicle 

V 

Solving equation (8-27) for the orbital rate gives 


(8-28) 


However^ for a vehicle on the surface of the earth the follow 
applies 


= '^<'3 

11m 


(8-29) 


where 


g = Gravitational acceleration on earth^s surface 


hence 


- r\ R,. 


(3-30) 


SubstltuLl.ng equation (8—30) into equation (8-28) yields 




(8-3/) 


Combining equation (8-31) with, those describing the angular 
command position, rate, acceleration, and rate of change of ac- 
celeration profiles give these in terms of the orbital altitude 
and time only with all other terms being constants* 

For ease In initializing the computer simulation for the 
Coincident Gimbal to perform earth tracking it is desirable to 
have orbital angle "0" as a function of telescope look angle 
This is accomplished by solving equation (8-*22) for ''9” as a 
function of 4> which results in 




+ C/J WV’} 

8.2 Comparison Between Actual Earth Trackine Profiles and 


the Hyperbolic Approximations - Approximate earth tracking pro 
files were employed in the evaluation of earth tracking perform- 
ance of the lOG and the Coincident Gimbal system for the suspension 
parameters specified by Dbrnier as described in sections 4i 2 and 
7.2 respectively. In. order to show that the results obtained 
using these profiles are valid, a comparison between the actual 
prpviles derived is section 8,1 and the hyperbolic approximations 
des bribed in s ectibh 4 . 2 vrill be p er formed in what f o Hows, 


c-i 



In order to perform tills comparison, the orbital altitude 
for which the approximations are -valid -must be determined* In 
order to accomplish this a relationship that relates peak track- 
ing rate to altitude must be determined. This can be done by 
noting that peak angular tracking rate occurs for 0=0. Hence 
from equation (8-19) the following applies 

( 8 - 33 ) 

Substituting equation (8-31) into equation (8-33) and re- 
arranging terms yields the following relationship 

- ■ 

=0 (8-39) 

yWx 

Solving equation (8-34) for the peak tracking rate of 1.748xl0' 
rad/sec when the approximate hyperbolic earth tracking profiles are 
employed yields an altitude of 438 KM. 

It is now necessary to determine the peak angular accelera- 
tion and the look angles for which they occur for the approximate 
profiles. In order to facilitate this determination the profiles 


given in section 4.2 are i<n:itten in functional form. 

- K ‘iar? ^ I (8 -3s) 

= “ Krt (8“36) 

Differentiating equation (8-36) twice yields the follovvring 

w ( l ) = (8-37) 

loii)- (K, roli) — seck^lk, ( 0 - 39 ) 

in order to determlrxe the peak acceleration point equation 
(8-38) is set to aero which results in the following 

Ca-39) 

■ where ' ■ ■ 


Substituting the exponential fo3nu for the hyperbolic functions 
in equation (8-39) yields the following 

-J -o (8-41) 

Solving equation (8-40) for e^ gives 

e’^ = (VT *-tfr)/z~ /,932 (8-42) 


^ - O. £, 585 " 


( 8 - 43 ) 


Using the values for and a defined in section 4.2 gives 
372-t 


bS 




( 8 - 44 ) 


or the time at which peak acceleration occurs is 

•i = 37?. 7. sec 


( 8-4 £ t ) 


Putting equation C8-37) in exponential form and substituting 
equation C8'*42) yields 




L732- 

0.93X 


I 

—y 


0.7i7SK<=^- ( 8 - 46 ) 


Using the values defined in section 4.2 for K and a gives the 
following peak acceleration 

~ 7.07)f/0-'’ f.ae//s«c^ (8-47) 

The telescope look angle for which this ocuiirs can he deter- 
mined by substituting the appropriate values into equation (8-35) 
which yields 

4' = aisy? rjc/ = SZS-a ^ 

The peak acceleration and the telescope look angle at which 
this peak occurs for the actual earth traclclng profiles are 1.856x10 
rad/sec^ and 31.1 degree respectively. These values compare favor- 
ably vrith those obtained using the approximate profiles thus valid- 
ating the results obtained using the approximations. 


8.3 Approximate Relationsliip Between Peak Angular Accelera- 
tion and Orbital Altitude — Using the relationships given in equa- 
tions (8-23) and (8—26) the peak angular acceleration as a function 
of orbital altitqde can be computed. A series of these computa- 
tions are given in table 8-1. 

The values shoxTO in table 8-1 are plotted in figure 8-3. 
Examination of this figure Indicates that the peak acceleration 
can be written as 

( 8 - 49 ) 

where b is the slope of the line shown in figure 8-3 XiJhich 
is -2.099. Using a least square fit to estimate the constant 
K yields the following 

a 

== fi Ce^so) 

where h is iii miles 

Equation (8-50) can also be x^ritten as 


where h. is in Km. These expressions can certainly be used with 
good precision (error less than 3 percent) over the range of their 
derivation and apply quite well to aiticud^s o£ 1000 miles, yield- 
ing r^.sults that are within 15 percent of the actuaJ. values. 

8»4 Approximate Relationships Between Peak Rate of Change 
of Acceleration and Orbital Altitude ^ Examination o£ equatiqh ; 
(8^24) indicates that peak rate of change of angular acceleration 
occurs when 0 equals zero^ Setting 6 equal to zero and substituting 
equation (8-28). yields 



Table 8-1. Peak Angular Acceleration as a Function 
of Orbital Altitude 


ORBITA 

L ALTITUDE 

PEAK ANGULAR ACCELERATION 

hii.es 

KILOMETERS 

RAD/SEC^ 

100 

160.9 

1.494x10’"^ 

150 

241.4 

6.478x10*"^ 
. “4 

200 

321.9 

3.557x10 

250 

402.3 

2.224x10 ^ 

300 

482.8 

1.508x10 ^ 

350 

563.3 

1.082x10 ^ 
. -4 

400 

643.7 

3.102x10 

-5 

450 

724.2 

6.258x10 
, ^ , -5 

500 

804,7 

. 

4.957x10 ■ 

— 










Using the above expression the peak rate of change of acceler- 
ation is given in table 8-2, 

A plot of the values given in table 8-2 is shovin in figure 
8-4. Exaiaination of this indicates that the peak rate of change 
of acceleration can be written in the same form indicated in 
section 8.3 for the peak acceleration. Following the same pro- 
cedure as outlined above, the peak rate of change of accaleratibh 
can be written as 

where h is in Km 

The approximate expressions derived in sections 8.3 and 8,4 
wil] be used in the explanation of the earth tracking performance 
of the Coincident Gimbal system. 


/. 832»/0 


fa-fi-z) 


ORBITAL ALTITUDE 
KM 


PEAK RATE OF CHAMGE OF 
ACCELERATION RAD/SEC^ 


160 

2A0 

320 

400 


2.245x10“^ 

6.415x10"^ 

2.611x10''^ 

1.29x10 

7.209x10"^ 


480 











PEAK AGCELERATION TIME RATE OF CHANGE (rad/sac 
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COINCIDENT GlilBAL SYSTEM EARTH POINT TRACKING PERPORMANCE 


In this section the earth tracking performance of the Co- 
incident Gimbal system is described. The torque required to per- 
form earth point tracking is estimated analytically and compared 
to the values obtained from the time simulation. An explanation 
of the large difference between the estimated and actual torques 
required is presented and verified via subsequent computer simu- 
lation. The payload mass and inertia that can be accommodated 
by the IPS in the earth point tracking mode as a function of 
orbital altitude is then established under the constraint the 
gimbal torque required will not exceed 80 percent of the torque 
capability of the gimbal torquers, or 16 n-ra. The tracking ac- 
curacy, pedestal rotations and translations with the accompanying 
isolator elongations, is also specified. All of the earth point 
tracking studies described in this section were performed for the 
recommended suspension parameters given in table 7-1. 

9.1 Estimated Torque Required to Perform Earth Point Track- 
ing - In order to determine the telescope mass and inertia that 
could be accommodated by the IPS in the earth point tracking mode 
under any gimbal torque constraint, the torque required to perform 
earth point tracking should be estimated analytically in order to 
give a reasonable starting point in this determination., 


This can be done by using the approximate relationship for 
peak angular acceleration shox'm in equation (8-50) for the par- 
ticular orbital altitude of interest. An examination of orbital 
altitude requirements of experiments presently anticipated for 
integration with the IPS show that they vary between 160 and 480 
Km, Choosing the 160 Km altitude, the approximate peak angular 
acceleration that will be encountered by the experiment is given 
by 


u 


64.06 


peak 


(160) 


2.099 


=1.514x10 ^ rad/sec^ 


(9-1) 


The peak torque required to perform earth point tracking 
should be between the torques required to achieve the peak angular 
acceleration in equation (9-1) when the telescope rotates about 
its center of mass to when it rotates about the hinge point. Using 
the mass properties for the nominal telescope given in . table 2—1 
the telescope inertia about its CM assuming that the earth track- 
ing profiles will be applied about the y axis only is 


(9-2) 


J_ = 2.648x10^ Kg-m^ 

3y 

The y axis moment of inertia of the telescope about the gimbal 
hinge point is given by 

J»y = 2.648x10^ + (2.683x10^) (1.689)^ = 1.03x10^ Kg-m^ (9-3) 

Hence, the torque required to perform earth point tracking 
from a 160 Km orbit should be between the minimum and maximum 
values listed below 

T . = (2.648x10^) (1.514x10"^) = 4.009 n-m (9-4) 

T = (1.03x10^) (1.514x10”^) = 15.59 n-m (9-5) 

max 

IPS earth point tracking performance for the nominal telescope 
from a 160 Km orbit is shoim in figures 9-1 through 9-21. Examina- 
tion of y axis control torque shoxiui in figure (9-20) indicates that 
the peak gimbal torque required for this earth tracking maneuver 
was 47 n-m, a factor of three larger than the peak required torque 
estimated above. 

A possible explanation of this phenomenon is that the inte- 
grator in the control torque formulation is charging to too large 
a value hence causing the large overshoot in y axis control torque. 

An examination of figure 9-20 does not seem to bear this out since 
an overcharged integrator should show up as a torque spike on the 
time scale of figure 9-20 since the pointing control loop bandwidth 
is two Hz and the integrator would discharge rapidly. However, 
in order to substantiate that the integrator in the y axis control 
torque formulation is not causing the observed large y axis control 
torques, IPS time response was determined for the same conditions 
as above only with the integral gain in the y axis controller set 
to zero. Figures 9-22 through 9-42 summarize these results. Com- 
parison of these results with those obtained above (i.e,, figures 
9-1 through 9-21) indicate that the only differences in response 
were in the y axis tracking position and rate errors (figures 9-23 
and 9-26 respectively) which, as expected, was appreciably larger 
than those incurred when the y axis integrator gain was at its 
nominal value. All other system parameters including the y axis 
control torque, exhibited virtually identical responses as those 
obtained when the y axis integrator gain was at its nominal value, 
thus substantiating that the integrator in the y axis controller 
was not the reason for the large y axis control torque. 


Another possible explanation for the above phenomenon is that 
the IPS pedestal is applying hinge forces to the telescope due to 
its dynamic behavior xfhich, in turn, cause disturbance moments 
that must be counteracted by the telescope control system. In 
order to determine the validity of this contention, the telescope 
hinge point and pedestal CM were made coincident with the IPS sus- 
pension center of elasticity and the computer run for landmark 
tracking from 160 Kin orbit was repeated. These revisions should 
minimize the dynamic hinge forces applied to the telescope by the 
pedestal since the rotation about the suspension center of elas- 
ticity Is also consistent with rotation about the pedestal CM 
and telescope hinge point , thus greatly reducing hinge point 
linear accelerations and the accompanying telescope disturbances. 

The IPS time response for these conditions are shown ip figures 
9-43 through 9-51. Examination of these results indicate that the 
torque required to earth point track was 15.9 n-m (see figure 9-51) 
which is in line with the analytical projections. The y axis con- 
trol torque peak of 23.8 n-m that occurs at 3.7 seconds into the 
run is primarily due to a transient condition In the Initialization 
of the run since not all states are set at their proper values 
corresponding to when the earth point tracking profiles are 
picked up. Therefore, the second peak which occurs at 28.5 sec- 
onds into the run is indicative of the torque required to earth 
point track. These results indeed verify the above contention 
that the pedestal dynamic motions are resulting in hinge forces 
which give rise to telescope disturbance torques which must be 
overcome by the control system. Hence, from the earth point 
tracking viewpoint, it is desirable that the telescope hinge 
point and pedestal CM be as close to the Suspensibh center of 
elasticity as possible, within the constraints of practical de- 
sign considerations. ■ 

However, previous studies (i.e., -'Evaluation of Alternate 
Telescope Pointing Schemes" NASA Contract No . NAS8-30889) indi- 
cate that making the telescope hinge point and the pedestal CM 
coincident with the center of eiasticity of the IPS suspenslpn 
has a detrimental effect on pointing performahce for disturbances 
applied in the X,y plane. Figure 9-52 is a reproduction of fig- 
ure 3-4 from Volume I of the final report titled, "Evaluation 
of Alternate Telestope Pointing Schemes,'' issued May 1975. Ex- 
amination of the curve for which r. = R = R_,n = 0^ which 

1 e40 340 

makes the pedestal CM and telescope hinge point coincident with 
the suspension center of elasticity, Indicates, that if the linear 
suspension stiffness exceeds approximately 4o6 n/rn^ 1 arc-second 
pointing stability trill be exceeded in the presence of crew motion 





l:u 


disturbances for a zero degree telescope look angle. The reason 
for this phenomenon is that the telescope hinge point does not 
translate under the Influence of the applied disturbance and the 
control torques applied to tne pedestal by the pointing control 
loop, which reduces the translation acceleration disturbance coup- 
ling into the telescope. In fact, with the recommended linear sus- 
pension stiffness of 2000 n/m to reduce pedestal translation and 
rotation incurred during landmark tracking and improve stellar 
pointing performance, the resultant pointing error for zero tele- 
scope look angle for a crew motion disturbance along the orbiter 
y aais obtained from figure 1 is approximately 5.3 arc-second. 

If pointing stability requirements are to be achieved, the sus- 
pension linear stiffness in the x and y directions would have to 
be reduced by an order of magnitude (i,e., 200 n/m). Although 
this reduction in linear stiffness in the x,y direction would 
tend to increase the pedestal translations incurred during earth 
point tracking, making the gimbal hinge point, pedestal GM, and 
suspension center of elasticity to coincide reduces the pedestal 
translations and rotations incurred. Hence, a trade should be 
conducted which determines the pedestal translations, rotations, 
and isolator elongations for different separations betoeen gimbal 
hinge point, pedestal CM, and suspension center of elasticity 
once the suspension parameters have been set to yield satisfactory 
stellar pointing performance for each particular set of geometrical 
parameters being considered. This type of study would indicate 
if there is a geometrical configuration coupled with a set of 
suspension parameters that would yield better overall (i.e. > 
stellar pointing and earth point tracking) IPS performance than 
those realized with the specified geometry and recommended sus- 
pension parameters , 

9.2 Payload Mass Characteristics That Can Be Accommodated 
By the IPS In the Earth Point Tracking Mode for Various Orbital 
Altitudes - In the section that follows, the payload mass charac- 
teristics that could -be accommodated by the Coincident Gimbal 
system, in the ^rth point tracking mode is described. In this 
determination the nominal payload mass and Inertia values given 
in table 2-1 were varied about all three axes by the same ratio 
while keeping the distance from the hinge point to the telescope 
(i. e. , rp constantv In addition, the control gains Were ad- 
justed in each case to maintain a pointing control loop bandwidth 
of two Hz* The gimbal torque was constrained not to exceed 80 
percent of its full torque capability of 2:0 n-ra which set the 
maximum allbwabie gimbal torque at 16 n-m, during the deter^na- 
tion of the payload mass and inertia characteristics that could 
be accommodated by the IPS during earth point tracking; 
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For the conditions enumerated above the telescope mass and 
inertia characteristics that could be accoramodated by the IPS 
during earth point tracking as a function of orbital altitude a.nd 
the peak tracking error incurred are given in table 9-1* 

The values given in table 9-1 are plotted in figure 9-53. 
Examination of this figure indicates that the inertia that can be 
accommodated by the IPS during earth point tracking without ex- 
ceeding a required gdmbal torque of 16 n-iu can be writtep. in the 
form 

(9-6) 

where b is the slope of the line shorn in figure 9-53/ Measuring 
this slope gives a value of 

h ^ Z.151 y (9-7) 


This value is trithin 2,7 percent of the negative slope of 2«099 
shown in equation (8-50) which gives the approximate relatiDnship 
between the peak acceleration and orbital altitude. Hence it is 
seen that the telescope inertia that could be accoimnodated by the 
IPS is a funetioh of the peak angular acceleration required to 
earth point track. Or another way of looking at this result is 
by realizing that the peak torque required is proper tipnal to the 
telescope inertia times the peak acceleration or 

^peak ^ peak 

Substituting equation (8-50) into equation, (9-8) gives 

,,23:,.6K^J,;, , , 

^peak " ^2.099 

Solving equation (9-9) for the inertia keeping the peak torque 
constant at 16 n-m gives the following 

: ^ ' ■ . 

which is within 2.7 percent of the results obtained above and 
defined by equations (9-6) and (9-7) . 


Table 9-*l. Telescope Tracking Error Mass and Inertia Characteristics 
That Can Be Accommodated by the IPS During Earth Point 
Tracking 


ORBITAL 

ALTITUDE 

(Km) 

TELESCOPE 

MASS 

(Kg) 

TELE8C0PE 

IHERTIA 

Kg-ra^ 

TELESCOPE Y AXIS PEAK 
TRACKING ERROR 
(rad) 

TELESCOPE Y AXIS PEAK 
TRACKING ERROR 
(ARC-SECOND) 

160 

1.309x10^ 

1.292x10^ 

9.93x10"^ 

20.47 

■ ■ ■ ■ 

240 

2.937x10^ 

2.898x10^ 

2.55x10“^ 

5.258 

320 

5.698x10^ 

5.622x10^ 

1.11x10“^ 

2.289 


8.987x10^ 

8.867x10^ 

5.76x10"® . 


480 

1,422x10^ 

1,402x10^ 

3.31x10"® 











9.3 Peak Tracking Error Incurred During Earth Point Tracking - 
Figure 9 t' 54 is a plot of the va3,ues given in tsbie 9-i for tele- 
scope tracking error vs orbital altitude incurred during earth 
point tracking. Examination of figure 9-54 indicates that the 
tracking error can be expressed in the same functional exponential 
form as the tielescope inertia discussed in section 9.2. The 
exponent which is equal to the slope of the line in figure 9-54 

is measured to be 3.092. This value compares favorably with the 
3,136 exponent in the expression for the rate of change of accel- 
eration given in equation (8-52) being within 1.4 percent of it. 
Hence j it is seen that the telescope tracking error incurred 
depends on the rate of change of angular acceleration which is 
directly proportional to the rate of change of torque rather 
than the peak tbrque required to perform earth point tracking. 

This is not surprising since there is Integral control and the 
effectiveness of the integral control in reducing telescope 
tracking error depends on how fast the torque requirements change. 
This implies that if the pointing control loop bandwidth were 
kept constant, regardless of telescope inertia variations, then 
the peak tracking error incurred would similarly remain constant. 
This would be exactly true if it weren't for the effect pf the 
disturbance torque applied to the telescope due to pedestal dy- 
namics which is appreciable as described in section 9.1. Table 
9-2 lists the telescope tracking error incurred for a 160 Km 
orbit as a function of mass and inertia variations for the re- 
commended suspension parameters and a 2 Hz pointing control loop 
bandwidth. Examination of table 9— 2 Indicates that telescope 
tracking error indeed does not follow the variation in telescope 
mass and inertia. The tracking error incurred for the nominal 
telescope mass and inertia (i.e. , 2.683xlo3 Kg, 2.648xlo3 Kg-ra^) 
and for one-fifth those values (i.e., 536.6 Kg; 529.6 Kg-m^) are 
essentially the same. This essentially substantiates the con- 
tention that the telescope tracking error incurred is essentially 
independent of telescope mass and inertia or required torque for 
a given orbital altitude as long as the pointing control loop 
bandwidth is kept constant. The relatively large pointing error 
incurred at prte-i^half nominal telescope mass and inertia is con- 
sidered to be due to a peculiar interaction with the pedestal dy- 
namics.^.. ' 

9.4 Fedestdl Hotatiohs Translations and Isolator Elongations 
Incurred During Earth Point Tracking - Table 9-2 gives the pedestal 
peak y axis rotations hnd x axis CM translatlpns incurred during 
earth point traclcing for the cases discussed above. 


Table 9“2, Pealt Pedestal Rotations and CM Translations Incurred 
Curing Eartb Point Tracking 


PEDESTAL PEAK Y AXIS 
ROTATION (rad) 

PEDESTAL PEAK Y AXIS 
ROTATION (deg) 

PEDESTAL PEAK X AXIS 
CM TRANSLATION (mm). 

PEDESTAL PEAK X AXIS 
CM TRANSLATION (in) : 

.■■■":0;^ll4l 

' :vV 

65.76 

2.589 

o^iiss;'. 

:> 6y52f; 

?/ 65.94 . 

2.596 

- bvine;. 

, /r ^6^394" 

65.89 

:: 2.594 

• 0;.1099 


66.45 

2.616 

d;i095 , 

^:6,274;; . 

:i;v 67,76 

2.668 


3 c:,-3 




ri‘ 













Examination of table 9-2 indicates that the pedestal rota- 
tions and corresponding CM translations are essentially constant 
for the cases investigated. Since the peak torque required in 
each case was kept constant at 16 n^-m, this implies that the 
pedestal trans lations and rota.tions depend to first order on the 
peak torque required to earth point track, and does not depend 
on the rate of change of the required torque (i.e,, its 
frequency characteristics) , This is the case since the rate 
of change of torque required to earth point track is still well 
below the bandwidth of the pointing control loop and isolator 
dynamics so that the angular rotations and CM translations of 
pedestal reach the same peak values although the points in time 
which these values ; occur change as the orbital, altitude changes . 

The isolator elongations associated with the pedestal rota- 
tions and translations given in table 9-2 is approximately 33 
ram Cl* 29 in) which is still too large to. be accommodated by the 
Coincident Gimbal suspension system design. Hence the telescope 
inertia that could be accoramodated during earth point tracking 
is limited by allowable isolator elongations rather than the 
gimbal torque limitstioh, and a deterrainat ion should be made of the 
telescope inertias that could be accommodated at various altitudes 
while yielding acceptable isolator elongations, 

9.5 Techniques for Improving Telescope Tracking Accuracy 
During Earth Point Tracking - It has been shown in section 9,3 
that the incurred telescope tracking error is primarily dependent 
on orbital altitude and essentially independent of telescope in- 
ertia as long as the loop bandwidth remains constant. Referring 
to figure 9-54 it is seen that a 1 arc-second peak tracking error 
during earth point tracking can only be maintained if the orbital 
altitude is greater than 420 Km for a 2 Hz pointing control loop 
bandwidth using angular rate and position commands. Since many 
earth pointing experiments are slated to fly at lower altitudes 
a method for improving tracking accuracy is desirable. One possi- 
ble way of improving earth point tracking accuracy is to add a 
torque comraani to the angular rate and position; cbramands already 
being implemented. One possibility that comes to mind is to add 
a torque command whose functional form will be the same as the 
orbital acceleration profile given in equation (8-23). However, 
as pointed out in section 9.1, the constant that would multiply: 
the orbital acceleration profile would have to be greater than the 
telescope inertia about the hinge point, A method of determining 
the constant that would multiply the acceleration profile would be 
to set that constant times the peak ac celerat ion at the time the 
peak torque occurs equal to the peak torque that is required to 
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perform earth trackdjig for a particular telescope at a particular 
altitude. It is clear that iJ this technique is to trork then the 
control torque reqviired to earth point track cannot deviate from 
the command torque which is functionally the same as the accelera- 
tion profile to any "large” extent* It should be noted that not 
onjy is the magnitude of the difference importaiit but the rate 
at which the torque divergence talces place is of critical im- 
portance since that determines the degree of effectlvity of the 
integral control. In order to determine just what the deviations 
in the actual required control torque are from the command torque 
which is characterized as described above, a 160 Km orbit with 
telescope and inertia mass characteristics that were one-half 
the nominal values described in table 2-1 and the recommended 
suspension parameters was chosen as a test case. Figures 9-55 
through 9-75 show system response for this case. 

Using equation (8-23) the torque command profile can be 
written as 



_ R h(R th) (2R th)sin[w t-0 ] 

f' 'S' • S S O' • O 

y i— — 

^ I A-n I \ « 


IR +(R + h) -2R (R+h)cos(uJ t-0 )] 
*■ e e e e ' o o' 



(9-11) 
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.'whuare 

0 - OrMtal angle with resp^^ to local vertical where earth 

point tracking coinmences . 

J„ = Constant that multiplies the angular acceleration profile 

For the 160 Km example cited the following numerical values 
for the terms in equation 9-11 apply 

' ' . ^3 ^ ' 

0) - 1.195x10 rad/sec 
o 

0 = 2,036x10 ^ rad /sec 

, o ^ ^ , 

R'-'-= 6378' Km' ■ " - ^ , 


n 





The point in time at which the peak control torque of -16.4 
n— m occurs is at 25.508 seconds into the run. However, the peak 
acceleration that occurs in the neighborhood of this time actually 
occurs at 29 seconds into the run. Hence, if the constant 





j is detertulned by requiring that the command torque and the actual 

I i required torque match at the point in time that the peak control 

torque occurs, the command torque xrould still be increasing until 
29 seconds is reached. Performing the actual computations the 
I angular acceleration at 25,08 seconds, the point at which the peak 

) control torque of -16.4 n-m occurs, is equal to 1,397x10“^ rad/ 

sec. The peak angular acceleration which occurs at 29 seconds 
j : is equal to 1.512x10“^, Hence if equation (9-11) were used as 

I ! the command torque profile the torque that would be commanded at 

29 seconds would be equal to 

1 \ -3 

^ T { = 16, 4 J = 17.75 n-m (9-12) 

29sec 1.397x10"^ 

1 ’ 

[. i The control torque is actually required at 29 seconds is 

-12.8 n-m or a difference of 4.95 n— m. This torque difference 
that must be made up by the control system is appreciable but 
perhaps even more importantly this torque difference occurs over 
four seconds of time. This makes the rate of divergence between 
the actual required torque and the command torque profile given 
j i by equation (9-11) appreciable which means that the integral 

' ■ control will not be very effective. The same behavior between 

peak acceleration and peak control torque is exhibited at the 
!j other altitudes investigated. Thus, although it is anticipated 

|j that using a torque command profile as described by equation V; 

(9-11) would reduce the telescope tracking error than those that j 

- - are presently incurred without it, it is doubtful that a tracking ■ 

I ^ accuracy of + 1 arc-second could be achieved. If such accuracies 

are desirable then the torque command profile would have to be I 

more complicated than that shoxro. in equation (9-11). It is I 

I therefore recommended that system performance with the command j 

profile shown in equation (9-11) and the form of the torque pro- ^ j 

files that would be required to yield + 1 arc-second tracking i 

errors should be the subject of future study. .1 
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10. INSTRUMENT POINTING SYSTEM (IPS) NOISE ANALYSIS 

This section will develop, via analytical techniques, an in- 
sight into the pointing errors as a function of control loop band- 
width for the IPS resulting from sensor and actuator (i^e. , gimbal 
torquers) noise. Various alternate sensing implementations are 
considered in these evaluations and Include the following: 

a. Separate Rate Plus Position Sensors 

b. Derived Rate 

c. Rate Gyro Control 

The torque range used by sensor noise for the sensor imple- 
mentations considered is also derived. These calculations will 
be used to determine the allowable sensor noise from the viewpoint 
of both pointing accuracy and the utilization of torquer range. 

In addition, optimal control loop bandwidths which result in 
minimum pointing error in the attitude estimate will be deter- 
mined as a function of gyro and position sensor noise levels. In 
this evaluation both white and colored noise for the gyros will 
be considered. 


10.1 Separate Rate Plus Position Sensor Implementation 


10.1.1 Pointing Error Due to Position Sensor Noise - A simpli- 
fied single axis linear representation of the telescope control loop 
using the inside-out gimbal configuration is shown in figure 10-1 
where 

Kp = position gain (n-m/rad) 

= rate gain (n-m/rad/sec) 

2 

J = telescope inertia (Kg-m ) 

N = sensor noise input 
T^ = control torque (n-m) 

The transfer function between sensor noise input N and angular out- 
put is given by 
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which can be ^jritten as 
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Equation (10-2) can be rewritten as 
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Expanding equation (10-5) results in 
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Factoring equation (10-6) gives 
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It can be shown that the mean squared vehicle pointing error 
is given by 
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where 


<i>i^(u)) = input noise power spectral density 
G^(t) = vehicle angular error 

Assuming a flat power spectral density (i.e., white noise) and sub^ 
stituting equation (10-7) into equation (10-8) gives 


-Hii 


OJ 


(tD+UJ^+j?w^) (w-£U^-jCw^) (10 9) 

where ^ is the amplitude of the flat power density spectrum. The 

integral shown in equation (10-9) can be evaluated via contour integra- 
tion using the contour shown in figure 10-2, 
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which can be written as 
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The rms pointing error Is given by 
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Assuming that 20 percent of the total pointing error will be 
budgeted for position noise the amplitude of the allowable posi-- 
tion noise power spectral density is given by 
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(10-15) 




where a damping ratio of was assumed. 


10.1.2 Angular Error Due to Rate Gyro Noise Assuming Separate 
Rate and Position Sensors - Using the block diagram shown in figure 10-3, 
the transfer function between position error and noise input is given 
by 
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Using the same techniques as outlined in section 10.1.1 
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Q = amplitude of flat rate gyro noise power spectral 


density 


/radV 

(sec, 

rad 

sec 


Assuming that 20 percent of the allowable 1 sec rms pointing 
stability error were alloted to rate gyro noise the amplitude of 
the rate gyro power spectral density is given by 
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K„ „ =5.323xl0"^^f f— 1 /radlsec 
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10.1.3 Angular Error Due to Torquer Noise - Using the block 
diagram shown in figure 10-4 the transfer function between angular 
error and torquer noise input is given by 
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Using the same techniques as outlined in section 10.1.1 the 
mean squared noise error is given by 
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where 


^ - flat power spectral density amplitude of gimbal 

(p-m) ^ 


torquers 


rad 

sec 


Assuming that 20 percent of the allowable 1 sec pointing 
error is budgeted for this noise source the following results 
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where a damping ratio of was assumed for the control loop* 


10.1.4 Torquer Range Used Due to Position Sensor Noi se - 
Using the block diagram shown in figure 10-1 the transfer function 
between sensor noise input and the resulting torque noise is 
given by ^ 

KpJs 
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Equation (10-24) can be rewritten as 
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Expanding equation (10-27) results in 
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which can be written as 
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Using the results shown above (equations (10-3) through (10-7) 
equation (10-29) can be written as 
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The mean square torque noise is given by 
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Examination of equation (10-31) indicates that if a completely flat power 
spectral density is assumed extending from Liie first Integral is 
unbounded as must be the case since at frequencies higher than the 


I 

"f 


10-6 


servo loop bandwidth the torque transfer function approaches Kp 

or Therefore a noise cutoff will be assumed which actually 

is the case since the sensors and actuators do have finite band- 
widths which were neglected in this simplified analytical model. 

This noise cutoff will only be assumed in the evaluation of the 
first integral of equation (10-31), The second integral will be 
evaluated from ^ as indicated. The second integral in equation 
(10-31) can be evaluated by contour integration using the path shown 
in figure 10-2, Therefore 
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where is the noise or sensor cutoff frequency. Simplifying 
equation (10-32) results in 
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Assuming that the sensor cutoff is five times the loop band- 


width then equation (10-33) can be written as 
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Assuming a damping ratio of equation (10-34) can be written 


as 
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Equation (10-35) can also be written as 


T^(t)=7.617xlO^J^f\ 
o n p.s. 


(10-36) 
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Substituting equation (10-15) into equation (10-36) yields 


T^(t)=5. 135x10 -^J f 
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The star sensor noise equivalent angle as a function of allowable 
torque noise can be obtained from equation (10-36) by assuming that 
the sensor has a bandpass that is five times greater than the 
corresponding vehicle loop bandwidth which can be expressed by 
the following relationship 
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where 


- sensor noise equivalent angle 

i (u) « sensor noise power spectral density 
ss 

tlgCco) = sensor transfer function 


Substituting equation (10-36) into equation (10-39), rearranging, simpli- 
fying, and taking the square root results in 
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The position sensor noise equivalent angle can also be deter- 
mined using equation (I*" -15), This results in the following 
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It should be noted that the position sensor noise equivalent 
angle was computed over a frequency range equal to five times the 
IPS pointing control loop bandwidth while using the value of power 
spectral density amplitude determined in equation (10-15). This was 
done in order to more realistically compare the impact resulting 
from pointing stability and torquer range utilization* 
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Equation (10‘-33) indicates that the sensor noise equivalent 
angle must be constant for any control loop bandwidth. It should 
be noted, however, that if the position sensor is used in an up- 
date mode the noise equivalent angle shown in equation (10-33) would 
only have to be met over the position update loop bandwidth* This 
can result in a large alleviation of the position sensor noise 
requirements as the position update loop bandwidth is decreased. 


The value of loop bandwidth for which the position sensor 
noise equivalent angle is equal from both allowable torquer range 
utilization and pointing stability considerations is given by 
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10,1,5 Torquer Range Used Due to Rate Gyro Noise Assuming 
Separate Rate and Position Sensors - Using the block diagram 
shown in figure 10-3 the transfer function between rate noise 

input and torque output is given by 
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Using the same technique as outlined In section 10,1,4 the 
torque noise resulting from rate gyro noise is given by 
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Assuming that the sensor bandwidth is five times the IPS point- 
ing control loop bandwidth equation (10-45) can be written as 
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Substituting equation (10-19) Into equation (10-47) results in 

T (t)i =1.435xlO“Sf^ (n-tn) (10-48) 

o * rms n 


■^n 

Ji 


The rate gyro noise equivalent rate error assuming that the 
rate gyro loop bandwidth is five times larger than the IPS loop 
bandwidth is given by 
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Substituting equation (10-47) into equation (10-49) gives 

.4 




0.1276 


R Jf 


T (t)| 


rad 
rms sec 


2,631x10' , . I sec 

Jf o^ ^'rms sec 

n 


(10-50) 


The rate gyro noise equivalent rate error without considera- 
tion of resultant torquer range utilized can be obtained from 
equation (10-19) as 
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where the integration was performed over five times the IPS pointing 
control loop bandwidth. 

The value of loop bandwidth for which the rate noise equivalent 
error is the same from both torque noise and pointing accuracy view- 
points is given by 
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Figures 10-5 through 10-10 summarize the results obtained in 
10,1,1 through 10*1,5, 

10,1,6 Computation of Optimum Pointing Control Loop Bandwidth 
for Separate Rate Plus Position Sensor Implementation - Using the 
results obtained in sections 10.1.1 through 10.1.3 the total mean 
squared error incurred from positions rate, and torquer noise can 
be written as 
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Differentiating the above expression with respect to the IPS 
loop bandwidth f^ and setting the results equal to zero gives 
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Equation (10-54) can also be written as 
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Solving the above equation for results in 
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Equation (10-56) expresses the IPS loop bandwidth that would 
result in minimum pointing error in the presence of position, 
rate, and torquer noise for the separate rate plus position 
implementation as a function of these noise levels. 

10.2 Rate Gyro Hold Sensor Implementation 


10.2.1 Pointing Error Due to Rate Gyro and Position Noiae - 
The simplified block diagram shown in figure (10-11) represents the 
rate gyro hold sensor Implementation. The transfer function 
between rate gyro noise input and IPS angular error is given 

fay 
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The power spectral density of the resultant pointing error 
can be i^nritten as 
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(10-59) 


Examination of equation (10-59) indicates that there is a double 
pole at the origin implying that an infinitely large mean squared 
angular pointing error would result if the rate gyro hold mode 
were implemented as shown in figure 10-11 « This result is not 
suprising since there is power around dc for the flat power 
spectral density of gyro rate noise that was assumed and hence 
the system will drift as time goes on resulting in ever increasing 
angular pointing error. It is therefore necessary that a gyro 
update scheme be employed in order to eliminate the effects of 
gyro drift thus keeping the angular pointing error finite. An 
analog update loop that can be employed for this function is 
shown in figure 10-12. 
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The position estimate power spectral density can be written 
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as 


Th^ mean squared error of the position estimate can be written 




0^(t)=1TtU K + — 
o u p.s. ai 

^ u 


(10-61) 


In order to find the optimum update frequency equation (10-61) 
can be differentiated with respect to and the results set 

equal to zero. This results in 

(10-62) 



or 
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,K 

V p.s. 
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The optimum mean squared error can be obtained by substituting 
equation (10-62) into equation (10-61) resulting in 
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0^(t)=2TT(K _ \ g 

o p.s. iv« G ^ 


(10-64) 
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and the xms error in the position estimate is given by 

.IM 




rms 
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(10-65) 


If the position sensor is truely updating the rate gyro then 

Therefore the mean squared error given in equations (10-61) 

and (10—64) is also the mean squared position error of the vehicle 
due to what can be thought of as an equivalent position sensor. 
However, the rate gyro still feeds the signal needed for rate 
damping into the system and to furnish the high frequency atti- 
tude information required for vehicle control. In order to ob- 
tain the gyro high frequency (i.e., a)^<a)<tiJ^) error contribution, 

equation (10-59) will be integrated from However, contour 

integration cannot be immediately applied on equation (10-59) since 
there is a double pole located at the origin. Therefore equation 
(10-59) is written in the following form 
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Upon evaluating constants A, B, C, and D the following results 
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Substituting equations (10^67) through (10-71) into equation (10’=-66) 
following results 
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Examnation of the first term In the brackets of equation (10-72) 
indicates that an infinitely large angular error will result when 
that term is integrated between This is to be expected since 

as indicated above it reflects the fact that the gyro will cause 
a vehicle drift since there is energy in the flat rate power spec- 
tral density assumed for the rate gyro around dc. However, the 
gyro drift error is taken care of by the gyro update loop as out- 
lined above. The high frequency noise contribution of the rate 
gyro can be obtained by integration of the second two terms in 
equation (10-72) over all of frequency. Performing the evaluation 
results in 
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Therefore the total vehicle mean squared pointing error due to 
sensors can be written as the sum of the low frequency errors 
given in equations (10-61) and (10-64) and the high frequency error 
given in equation (10-73) which results in 
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Substituting the optimum update frequency value given in equation 
(10-62) results in 
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The overall mean squared pointing error due to both sensors 
and actuators is given by 
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Examination of equation (10-76) indicates that there is no clear 
optimum loop bandwidth, for minimum pointing error^ at which the 
rate gyro hold loop should operate* The noise error incurred will 
decrease with increasing loop bandwidth and will be asymptotic to 


asump 

These results have been verified on various LST simulations where 
rate gyro hold was employed. 

Assuming that 20 percent of the allowable 1 sec pointing error 
would be budgeted for high frequency rate gyro noise the following 
results 



„_.4.236xl0-l\ (10-78) 

sac 

1 

where a damping ratio of was assumed. 

42 

Assuming that the position sensor and torquer noise power 
spectral density are those given in equations (10-15) and (10-19), 
respectively, the resultant mean squared pointing error for the 
rate gyro hold sensor implementation is given by 

1/2 2 

6^(t)=2fr(5.323xl0“^^x4. 235x10“^^) +2[ .2(4.85x10”®) 1 -1.132xl0"^^ rad^ (10-79) 

0o<t)|j.^tal"^*^^^^^'^~^ rad=.6936 sec (10-80) 

rms 

It is therefore seen that the total rms IPS pointing error 
will be met with the numbers budgeted for the various error 
sources . 


10.2.2 Torque Noise Due to Rate Gyro Noise for Rate Gyro 
Hold Sensor Implementation - Using the block diagram shown in 
figure (10-11) the transfer function between gyro noise input and 
the resultant torque noise is given by 
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Using the same techni^iues as outlined in sections lOal.4 and 
10.1.5 the mean squared torque noise resulting from rate gyro 
noise is given by 

»2 3„ 




(10-83)- 


Assuming a rate gyro cutoff frequency of five times the IPS 
pointing control loop bandwidth the following results 
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l2(t)=40;Vj\_G * ^ K^.5.[1*4C^-16 c'^1 


(10-84) 


becomes 


Assuming a loop damping of -p; the mean squared noise torque 

i2 


T^(t)=17.7au^J^K„ _ 
o n tt.G. 


(10-85) 


The gyro mean squared rate error assuming a cutoff of five 
times the loop bandwidth is given by 




( 10 - 86 ) 


Substituting equation (10-85) into equation (10-86), simplifying 
and taking the square root results in 


T (t)j 

o -.1194 — Jj 
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f— ^ 

\aec / 
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The rms rate error allowable from pointing considerations 
only is given by 

al=1.829xl0“'^f (10-88) 

K n 


Therefore the loop bandwidth for which the allowable gyro 
mean squared rate error is equal from both pointing accuracy and 
torquer range utilization viewpoints is given by 
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f =80.81- 
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These expressions are essentially equal to those shown In 
section 10.1.5 and hence the plots of those expressions also 
apply here. 
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10.2.3 Pointing Accuracies and Update Frequencies Involved 
When Using the LOG 540 Gyro In a Rate Gyro Hold Mode - Using the 
noise data for the LDG 540 gyro obtained by Martin Marietta 
Aerospace, end presented at the Seventh Biennial Guidance Test 
Symposium >'t tiolloman Air Force Base In a paper titled ’’Testing 
Technology for Fine Pointing Systems” by R. L. Gates, the low 
frequency power spectral density of the LDG 540 can be approxi- 
mated as 

*^18 2 

v|2 2.211x10 rad nos 

j(|)(to)| 5 (10-90) 

sec 


where adjustments were made to the data presented In the paper to 
obtain the units employed In this report, and to account for the 
double sided power spectral densities employed throughout this 
study. Since the data Is presented In terms of an angular power 
spectral density referred to as table motion when it is controlled 
by the LDG 540 gyro operating in the attitude mode, the gyro update 
loop considered in section 10.2.1 can be revised in the following 
manner as shown in figure (10-13). 

It should be noted that this formulation favors the LDG 540 
gyro since the noise due to the electronics which would be present 
In the strapdown mode is not accoimted for In equation (10-90). How- 
ever, it is presently felt that the electronic noise contribution 
In the frequency range of Interest is not substantial enough to 
grossly change the results obtained. 


The transfer function between rate gyro noise input 
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angular error output is given by 


B(.). -jlip (10-91) 
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and 
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The power spectral density of the angular error in the posi- 
tion estimate due to rate gyro noise is given by 

— 1 ft 

|e(,,)|2. (10-93) 

iJj 
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The mean squared angular error in the position eatiiaate due 
to rate gyro noise is given by 

q 2, . I _ Tr(2.211xl0~^^) (10-94) 

o'- ■'<r.g7 jd 

u 

The total mean squared error in the position estimate due to 
both rate gyro and position noise is given by 
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The optimum update loop bandwidth resulting in minimum error 
In the position estimate is given by 


to = 
u 



(10-96) 


and the total mean squared error in the position estimate is given 
by 
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^ ^ =271(2.211x10 xK 

total p.a 
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The total pointing error due to both rate gyro and position 
sensor noise can be written as [i.e.» equation (76)] 
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Using the value for position sensor noise given in equation 
(10-15) and assuming an IPS loop bandwidth of 2 Hz with a damping 

ratio of ~ the following results 

0Q(t)|totai=9-^52xlo”^^ rad^ (10-99) 
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(10-100) 
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rad=.2005 sec 

rms 

This value of rms pointing error is almost totally due to 
torque generator noise and the contribution of the position sen- 
sor noise has been virtually eliminated by the use of the LDG 540 
in a rate gyro hold mode* The gyro update frequency is given by 

f =5.142x10“^ Hz (10-101) 

u 

or 

x^“l,945 sec (10-102) 

This calculation indicates the reduced sensitivity to posi- 
tion noise that can be obtained using a rate gyro in an update 
mode* However, it seems that the LDG 540 gyro is considerably 
better than that needed for the lOG* 

10,3 Derived Rate Sensor Implementation 

10*3.1 Angular Error Due to Position Molse for Derived Rate 
Sensor Implementation - A simplified block diagram of the IPS 
control loop using a derived rate sensor implementation is shown 

in figure 10-14* The transfer function between position noise input 
and the resultant angular error is given by 
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Using the same method as outlined above the mean squared 
angular position error can be written as 
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The ratio of the mean squared angular error due to position 
sensor noise using the derived rate sensor Implementation to that 
incurred using the separate rate plus position sensor implementa- 
tion is given by 


0^t) 


derived 

rate 
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=4^+1 


(10-106) 


0^(t)i 

o ‘ separate 
rate plus 
position 


For a system damping of -p^ this ratio is three or approxi- 

u 

mately 1*73 times the rms angular error due to position sensor 
noise results vfhen the derived rate sensor implementation is 
used over that which would be Incurred for the separate rate 
plus position sensor iraplementationi 


10*3*2 Torque Noise Due to Position Noise for the Derived 
Rate Sensor Implementation - Using the block diagram shown in 
figure 10-14 the transfer function between position noise luput 
and torque noise output is given by 
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Equation (10-108) can also be written as 
l-t- 8 C^- 16 C^* 


iHOu.) 




(4 ^ 2 ) 


where 


A=[(2-4i;'^) (l+8C^-16c'^)-4i:^lw^ 


B=[ 16 C^- 8 C^-l]to^ 


(10-107) 


(10-108) 


(10-109) 


(10-110) 

( 10 - 111 ) 


in-20 





Using the same technique as outlined in section 10. 1.4 the 
mean squared torque noise can be written as 

r ^ 

4 + j K (10-112) 


2 17 2 

r(t)-8j^?V 

o n 


4C 


n q 


p*s. 


Assuming the position sensor cutoff is five times the IPS 
loop bandwidth equation (10-112) can be written as 
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The ratio of mean square torque noise for the derived rate to 
the separate rate plus position sensor implementation is given by 
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Assuming a damping ratio of 


—^equation (10-114) results in 
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'derived 

rate 
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= 21.86 


(10-115) 


or approximately 4.7 times the rms torque noise will result in 
the derived rate sensor implementation as that which results for 
separate rate plus position sensor implementation. 

Assuming a damping ratio of ^=jthe mean squared torque for 
the derived rate sensor implementation can be written as 


T^(t)=173.3jVK ^ =1.697xloS^f\ (10-116) 

o np*Sfe np.s. 
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The allowable position power spectral noise amplitude assuming 
that 20 percent of the allowable 1 sec pointing error would be 
budgeted for position noise is 


K 

p.s. 


2.247x10"'^^ rad^ 

f rad 

n 

sec 


(10-117) 


The star sensor noise equivalent angle as a function of allow- 
able torquer range utilization can be obtained by assuming the 
star tracker bandwidth is five times that of the IPS pointing con- 
trol loop. Hence 

a^=207Tf K rad^ (10-118) 

R n p.s* 

Substituting equation (10-116) into equation (10-118) and taking 
the square root gives 

■5 T (t)| -1 T (t)| ^ 

rtnr -.«"3 o ' nras , ^ o *rras /in ^^n\ 

a =6.085x10 rad=l. 255x10 sec (10-119) 

^ jr jr 

n n 


The allowable star sensor noise equivalent angle from a point- 
ing stability viewpoint is given by 

o’=l. 188x10“^ rad=0.245 sec (10-120) 


The value of IPS loop bandwidth for which equal noise equiva- 
lent angle requirements result from both pointing stability and 
torquer utilization viewpoints is given by 


f =71.57 
n 



(10-121) 


10.3,3 Computation of Optimum Loop Bandwidth for Derived Rate 
Sensor Implementation - The total mean squared pointing error for 
the derived rate sensor implementation Is given by 
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Differentiating equation (10-122) with respect to setting 
the results equal to zero and solving for f^ gives 
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which is the bandwidth for which minimum pointing error results. 


10.3.4 Criteria for Choosing Between Derived Rate and Sepa- 
rate Rate Plus Position Sensor Implementations - The criteria for 
when one chooses between separate rate plus position and derived 
rate sensor implementations can be formulated by the following 
inequality 

10^(t)L„ (10-124) 
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Substituting equations (10--53) and (10^122) iut’^ equation (10"*12A) 
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Simplifying equation (10-125) results in 


<l6TT^f!K 


n p.s. 


(10-126) 


When the rate gyro noise power spectral density meets the 
above inequality separate rate plus position sensor implementa- 
tion should be employed. 


Another way of viewing equation (10-126) is to write it in the 
following form 

f_> (10-127) 


■■n- 4 ttIK 
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or separate rate plus position should be employed over derived 
rate when the IPS loop bandwidth meets the inequality given in 
equation (10-127)* 


10, 3. 5 Criteria for Cnooslng Between Rate Gyro Hold and 
Derived Rate Sensor Implementations - The criteria for choosing 
between derived rate and J:ate gyro hold sensor implementation 
can be formulated by the following inequality 
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Substituting equations (10-76) and (10-122) into equation (10-128) 
results in 
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After some algebraic manipulations equation (10-127 can be 
written as 
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If the position sensor noise power spectral density meets the 
inequality shown in equation (10-130), derived rate should be employed 
over rate gyro hold. 


Another way to viev; the inequality shown in equation (10-130) is 
the followi 2 
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f < fl2C^-bH-4C(8C^-t-l) 1 

““ 2n(4C^+l) 
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(10-131) 


p.s. 


If the IPS Loop bandwidth meets the inequality constraint 
shown in equation (10-131) the\ derived rate should be employed 
over the rate gyro hold sensor implementation. 


10*4 Summary and Recommendations - The following are some 
of the results during this study phase: 

a. The derived rate sensor implementation will give three 
times the rms pointing error as the separate rate plus posi- 
tion sensor implementation for the same value of position 
noise power spectral density for an IPS control loop damp- 

1 

ing ratio of • This is true regardless of IPS pointing 

control loop bandwidth. 

b. The derived rate sensor implementation will use approx- 
imately five times the torquer range as the separate rate 
plus position sensor implementation for the same value of 
position noise power spectral density for an IPS control 

loop damping ratio of ■■ ^ . This is true regardless of IPS 
pointing control loop bandwidth. 
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c* For the derived rate sensor implementation. If 2 n-m of 
torquer range is budgeted to be used by sensor noise, which 
would correspond to 10 percent of the gimbai torquer range 
for a 20 n-m torquer, the allowable sensor noise power spec- 
tral density would be governed by this requirement for IPS 
loop bandwidths of 2 Hz and above. For IPS loop bandwidths 
below 2 Hz the allowable sensor noise power spectral densi- 
ties would be governed by pointing stability requirements r. 

It is recommended that a study be performed to determine the 
effects of sensor t'^nd actuator noise on overall IPS performance 
using the system three body performance model described above 
which were performed on a simplified linear second order model 
of the IPS using white noise spectra. Actual sensor and activa- 
tor noise characteristics including the effects of sampling and 
quantization which would not be handled in an efficient manner 
analytically would be included, thus allowing an optimum choice 
of sensor implementation. 
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Figure 10-5. Noitie Torqu<‘ Due to Position Sensor Noise 
(Equation 10-38) 
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Figure 10- f>. Comparison of Position Sojisor Noise Equivalent Angle Koquircments 
from PoinLlng SLabiliLy and Tor<|ui* Ntilsf" Viewpoints 
(Equations 10-40, 10-41) 
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Figure 10-7. IPS Loop Bandwidth for liflilch Po.sition Sensor Noise Equivalent 
Angle Requirements are Equal from Both Pointing Stability and 
Torquer Noise Viewpoints (liquation 10-42) 
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Figure 10-9* Comparison of Rate Sensor Noise Equivalent Rate Error 
Requirements from Pointing St 'ability and Torque Noise 
Viewpoints (Equations 10-503 10-51) 
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Figure 10-10. IfS Loop Bandwidth for Which l^oine EahiVAlent Rate Err qr 
Requirements are Equal from Both Pointing Stability and 
Torquer Noise Viewpoints (Equation 10-52) 
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Simplified Block Diagram for Rate Cyro 
Hold Sensor Implementation 
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Figure 10-12. 
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11. SUMMABY AI® CONCLUSIOHS 

The major conclusions of the studies performed can be sum- 
marized as follows : 

a. IPS Stability and performance is not sensitive to sus- 
pension damping characteristics. Damping can be varied over 
mde tanges (i,e,, reduced by an order of magnitude from a 
nominal damping ratio of 0 . 1 ) without any substantital change 
in system pbintini performance. 

b. IPS stability and performance is not sensitive to pedes- 
tal mass and inertia variations with pointing performance 
remaining essentially the same for pedestal mass ahd iner^ 
tia variations of plus and minus an order of magnitude from 
nominal values. 

c. IPS pointing pe^rformance is a function of telescope mass 
and inertia characteristics , A pointing error of approximately 
2.6 arc-seconds was incurred for a pointing control loop band- 
width of 2 Hz in the presence of crew motion disturbances for 
the telescope mass characteristics listed beloxg’ 

J 33 J. = 50 Kg~m^ 

= 50 Kg-m^ 

=35 Kg-m^ 

^^^m2' = ^290^Kg^ 

If telescopes in this class are to be integrated with 
the IPS suspension system, modification would be required 
if + 1 afc-secdnd pbihtihg stability is to be achieved. 

d. A pointing error of 0.895 arc-second was incurred for 
a pointins ‘control loop bandwidth of ? Hz, in the presence 
of crew motion disturbances using the suspension system 
parameters specified by Dornier Systems for the Coincident 
Gimbal system configuration. This pointing error which is 

: due to one source only is quite large and it is doubtful 

if + 1 arc-second pointing stability could be realized 
once all the other contributing error sources (e.g. , 
sensor and actuator noise, gimbal nonlinearities , sampling 
and quantization, etc.) are taken into consideration. 


e. Isolator ^ial elongations of approximtely 4. 3 nmi (li 7 
in) result during earth point tracking from a 438 KH orbit 
for the suspension parameters specified by Dornier Systems 
for the Coincident Girabal configuration* Kies e elongations 
which become worse as the orbital altitude decreases, are 
too large to be accommodated by reasonable suspension sys- 
tem designs, 

f. It is recommended that the suspension systeni parameters 
for the Coincident Gimbal be modified fromi those specified 
by Dornier co the following: 

K = K = 2G00 n/m Ic « k = 202.5 n-m/rad 

; X y , x y ■ 

K - 500 n/m k “ 3.242x10^ n-ra/rad 

z . , _ z 

D = D =475.2 n-sec/m 

■ X ■ . :y . 

D = 118.8 n-sec/m 

The values for the suspension parameters listed above 
result in the following isolator design characteristics 
for the Coincident Simbal suspension system designs: 

K = 750 n/m 

D = 178.2 n-sec/m 

a = 19.47 deg 

R = 0.9 ra 

The modified suspension not only reduces isolator 
elongations incurred during earth point tracking but also 
improves pointing performance oyer that which resulted 
for the Dornier suspension Specifications. ■ 

g. If the gimbal wire torque nonlinearity is accurately 
described by the characteristic shown in section 5, the 
jump in that wire torque characteristic would have to be 
0,05 n-m or less for the nominal telescope and 0.01 n-m 
or less for the 290 Kg telescope discussed above in item 
c, if a pointing stability of + 1 arc-second is to be 
achieved. These values will result In limit cycle ampli- 
tudes of 0.2 arc-second which is twenty percent of the 
total allowable error budget. 


d = d = 48,11 n-m-sec 
X y , 

d = 769.8 n-m-sec 
z 


h. The utaxiiaum telescope inertia that could Ise accontmodated 
during earth point tracking assuming that not more than 80 
percent of the girabal torque capability of 20 n-m be used j 
is proportional to the peak angular acceleration required 

to earth point track at a given orbital altitude. However, 
the actual torque required to earth point track consider- 
ably exceeds, particularly at lower orbital altitudes, the 
peak angular acceleration times the telescope inertia about 
the gimbal hinge point due to the action of pedestal dynam- 
ics. This effect can be reduced as the gimbal hinge point, 
pedestal CM, and suspension center of elasticity are brought 
closer together. These results are for the recommended 
Coincident Gimbal suspension parameters and do not take 
into account the isolator elongations Incurred during earth 
point tracking but only considers the 16 n-m torque limi- 
tation. 

i. The telescope tracking errors incurred are appreciably 
exceeding + 1 arc-second for orbital altitudes below 420 
KM for a pointing control loop bandwidth of 2 Hz when only 
angular rate and position commands are employed. It is 
clear that if + 1 arc-second tracking accuracy is to be 
achieved during earth point tracking, a feed forward torque 
command would be required in addition to the angular rate 
and position commands presently being employed. 

j. Telescope tracking error is proportional to the rate 
of change of angular acceleration required to earth point 
track and is essentially independent of telescope inertia 
provided that the pointing contrcl loop bandwidth remains 
constant. 


12 . RECOMENDED FUTURE EFFORT 


The following tasks are recommended to continue and extend 
the results obtained in this study in order to better understand 
and aid ESA in the design of the IPS. 

a. Deteannine the functional form of the command torques 
profile that would be required to improve system tracking 
accuracy as a function of command torque profile complexity 
with a + 1 arc-second tracking accuracy for orbital altitudes 
between 160 and 480 KM being a design goal. 

b. Determine the effect of parameterizing the distance be*- 
tween the gimbal hinge point, pedestal CM, and suepenaion 
center of elasticity, on system pointing and slewing perform- 
ance. In particular, investigate the torque required, and 
isolator elongations incurred, to perform earth point track- 
ing as a function of these dii3tances taking into consideration 
the suspension modifiiiations required to maintain + 1 arc- 
second stability during stellar pointing as the distance 
betvireen telescope hinge point, pedestal CM, and suspension 
system center of elasticity, decrease. 

c. Determine the maximimi telescope mass and inertia that 
could be accommodated as a function of allowable isolator 
elongation for various orbital altitudes. 

d. Determine the effects of sensor and actuator noise on IPS 
pointing performance. Establish the allowable levels for 
these noise sources as a function of sensor implementation 

and specify recommended implementation update frequencies, etc. 
as a function of sensor and actuator noise characteristics. 

e. Determine the effects of sampling and quantization on IPS 

pointing performance. Establish the sampling frequency and 
quantisation levels that would result in satisfactory system 
performance without compensating for these effects. If these 
values result in excessive software utilization investigate 
use of compensatioiv to reduce sampling frequertcy while still 
maintaining desired system performance characteristics. The 
effects of sampling and quantization on the errors incurred 
due to sensor and actuator noise characteristics should also 
be determined • . 



I . f , DeterraiB.2 the effect of telescope flexibility on over- 

j all IPS stability and pointing performance. Establish a ^ 

meaningful technique by which telescope flexibility character- 
istics could be specified so that if adhered to by the pay- 
load designer would result in satisfactory system stability 
I and pointing performance once integrated with the IPS. 


i 

I 



